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I.  GENERAL  INTRODUCTION 

Inis  work  was  undertaken  over  a period  of  three  years,  during  which 
t.  Live  research  in  nonlinear  optics  was  increasingly  concerned  with  both 
phenomena  of  higher  order  than  nonlinearities  and  with  applica- 

tions of  nonlinear  optical  techi  Lques  to  spectroscopic  problems. 

This  final  report  discusses  progress  made  in  tb i growth  and  applica- 
tion of  the  chalcopyrite  crystals  CdGeAs2,  AgGaSe2  and  AgGaS2  to  the  gen- 
eration of  tunable  infra.-ed  radiation.  The  bulk  of  this  report  is  work 

I 

that  comprised  R.F.  Begley's;  Ph.D.  thesis.  However,  the  Appendices  include 

1 

work  of  Mr.  M.M.  Choy,  Dr.  D.  Chemla,  Dr.  R.L.  Herbst  and  Dr . S . Ciraci. 

] 

The  recent  tremendous  need  for  tunable  infrared  radiation  reaffirms  our 
early  expectations  that  the  chalcopyrite  crystals  will  make  an  important 
contribution  as  new  nonlinear  optical  materials. 

A.  NONLINEAR  OPTICAL  PROCESSES 

The  first  part  of  this  report  considers  the  development  of  high  peak 
power,  short  pulse,  tunable  infrared  sources  based  on  nonlinear  optical 
processes  in  semiconductor  materials.  The  primary  laser  source  considered 
is  the  atmospheric  pressure  TEA-C02  laser,  whose  high  efficiency,  large 
output  powers  and  relatively  small  operating  cost  make  it  an  ideal  basis 
for  high  power,  tunable  infrared  devices. 


The  principal  nonlinear  material  considered  is  CdGeAs2,  a chalcopyrite 
crystal  with  42  m symmetry.  The  linear  and  nonlinear  optics  properties  of 
this  material  have  been  examine.!  by  several  workers.1'11  Due  to  its  large 


birefringence,  CdGeAs^  can  be  phasematched  for  second,  third,  and  fourth 
harmonic  generation  within  a large  part  of  its  transparency  range  from 

2.4  to  18  nm.  In  addition,  this  crystal  can  be  phasematched  for  difference 
frequency  mixing  of  a CO  and  COg  laser  to  generate  wavelengths  between 

11.4  and  16.8  ^m.  When  pumped  with  5.3  nm,  CdGeAs^  looks  promising  as  a 
parametric  oscillator  source  covering  the  range  from  7 to  18  ^m.^  Moreover, 
this  crystal  has  the  highest  figure  of  merit  for  second  order  nonlinear 
interactions  of  all  materials,  with  the  exception  of  tellurium.  Figure  1^ 
illustrates  the  transparency  ranges  and  figures  of  merit  of  all  nonlinear  . 

crystals  currently  in  use.  Finally,  CdGeAs2  has  a large  burn  density  of 

, 2 k 

33  MW/cm  for  a 200  nsec  long  pulse  from  a TEA-C02  laser.  These  last  two 

facts  point  to  the  possibility  of  doing  very  efficient  nonlinear  interactions 

in  quite  small  crystals.  We  have,  for  example,  measured  second  harmonic 

conversion  efficiencies  of  2#  and  conversion  to  third  harmonic  of  lo"6  in  a 

2 mm  long  crystal  using  a 50  KW,  200  nsec  pulse  from'a  C02  laser. ^ If  we 

g 

extrapolate  to  a 1 nsec  pulse  with  10  W peak  power  and  a 1 cm  long  crystal, 
we  predict  theoretical  conversion  efficiencies  of  60$  for  second  harmonic 
and  30$  for  third  harmonic  generation. 

Such  large  conversion  efficiencies  indicate  that  both  three  and  four 
frequency  parametric  processes  can  be  well  above  threshold  in  just  a few 
millimeters  of  material.  The  three  frequency  parametric  process  can  be 
phasematched,  as  mentioned  above.  The  four  frequency  parametric  process, 
using  two  photon  .5  of  10.6  nm  as  a two  frequency  pump  to  generate  both  a 
longer  and  a shorter  wave,  is  phasematched  only  very  close  to  10.6  nm. 

The  calculated  parametric  gains  of  several  nonlinear  materials  with 
a length  of  1 cm  and  an  input  pump  intensity  of  1 MW/cm^  are  shown  in  Fig.  2.^ 
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For  CdGeAs  , a 5 $ gain  at  1 MW/cm^  implies  a gain  of  approximately  1.5 

2 

at  30  MW/cm  which  is  just  below  the  burn  density  for  200  nsec  long 

8 2 

pi'ises.  For  a 1 nsec  pulse  with  10  W/cm  we  reach  a gain  of  5 in  a 

2 

1 cm  long  crystal.  Burn  densities  of  1 GW/cm  reported  in  p-typ’ 

0 

germanium  with  a 1.7  nsec  pulse  raise  the  exciting  possibility  that 
burn  densities  in  CdGeAs^  may  be  about  the  same  with  1 nsec  pulses. 

For  a 1 GW/cm  input  intensity  super-radiant  gains  of  nearly  50  in  a 1 cm 
CdGeAs^.  oscillator  crystal  are  possible. 

In  order  to  fully  evaluate  the  potentials  of  this  material  for  effi- 
cient nonlinear  optics  re  have  constructed  two  C0Q  laser  systems.  The 
first  device  consists  of  a pin  discharge  atmospheric  pressure  laser  oscil- 
lator , which  is  modelocked  with  a germanium  acousto-optic  modulator.  The 
modelocked  pulse  train,  consisting  of  twenty  1 nsec  wide  pulses,  each 
separated  by  10  nsec,  is  passed  through  an  electro-optic  selector  arrange- 
ment which  picks  out  a single  1 nsec  pulse.  That  pulse  is  in  turn  ampli- 
fied  in  a three  pass  C0^  amplifier  to  reach  a peak  power  of  10  - 10 J W. 

The  second  laser  device  is  a high  pressure,  1 meter  long  C0o  laser 
oscillator  which  operaten  between  5 and  10  atmospheres.  At  these  pressures 
the  individual  rotational  lines  are  so  pressure  broadened  that  they  overlap 
significantly.  By  making  a three  mirror  cavity,  with  a diffraction  grating 
as  a tuning  element,  the  output  of  this  laser  can  be  continuously  tuned 
across  the  9 to  11  p range  of  the  CO^  molecule. 

Four  aspects  of  the  work  performed  at  Stanford  are  covered  in  this 
report.  First,  the  crystal  growth  of  CdGeAs0,  AgGaSe^,  and  AgGaS,  are 
described.  Second,  the  design  and  operation  of  the  1 nsec  system  is 
discussed  in  detail.  Third,  a theoretical  treatment  is  provided,  both 


I 


f 

t 


for  higher  order  nonlinear  optical  processes  in  CdGeASg,  and  for  the 
origins  of  its  large  third  order  nonlinear  susceptibility.^  Finally, 
a description  of  the  experimental  work  performed  to  date  with  this  source 
is  given.  Saturation  mechanisms,  particularly  free  carrier  absorption, 
are  discussed. 

At  this  point  one  might  ask  what  applications  require  both  high  peak 
powers  and  tunability  from  a source  operating  at  wavelengths  from  2 to 
20  ^m.  We  can  separate  the  applications  of  such  sources  into  three  cate- 
gories: spectroscopy,  chemical  kinetic  studies,  and  optical  pumping 

studies . 


Perhaps  the  most  active  area  of  laser  spectroscopy  at  present  is  the 


study  of  vibrational  rotational  spectra  of  molecules,  many  of  whose  lines 
lie  in  the  near  infrared.  Such  spectra  can  be  used,  for  example,  as 

indicators  of  the  composition  of  the  atmosphere  by  remote  air  pollution 

9 lo 

detection  schemes.  ’ A second,  and  perhaps  much  less  explored  area  is 
two  photon  absorption  spectroscopy. 

Most  optical  spectroscopy  of  solids,  for  example,  has  been  done  by 
reflectivity  techniques.  Since  the  first  experiments  on  two  photon  absorp 


tion  near  the  fundamental  absorption  edge  using  visible  and  ultraviolet 

12-17  . . 

sources,  xt  has  been  clear  that  this  technique  allows  one  to  probe  the 

properties  of  bulk  materials  as  well  as  of  the  surface  states.  The  principal 
advantage  of  this  technique  is  that  each  photon  lies  within  the  transparency 
range  of  the  material.  Tunable  infrared  sources  then  should  have  consider- 
able application  in  two  photon  spectroscopy  of  semiconductors  whose  band  gap 
energies  lie  in  the  near  infrared.  Moreover,  high  peak  power  sources  should 
allow  studies  of  nonlinear  interactions  not  accessible  by  conventional 
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sensitive  to  the  shapes  of  the  energy  bands,  for  example. 

19 

Applications  of  lasers  to  chemistry  have  been  reviewed  by  Moore. 
Vibrational  and  rotational  relaxation  processes,  as  well  as  energy  transfer 
mechanisms,  can  be  studied  with  tunable  infrared  sources  with  pulses  in  the 
nanosecond  and  picosecond  regimes.  The  question  of  inducing  chemical  reactions 
by  selectively  exciting  vibrational  levels  has  received  much  attention, 
particularly  with  regard  to  isotope  separation. 

Finally,  optical  pumping  of  high  pressure  gases  appears  to  be  an 

important  technique  for  the  generation  of  wavelengths  across  both  the  near 

20-27 

and  far  infrared  with  pulse  widths  of  1 nsec  and  shorter.  With  high 

peak  power  tunable  sources  one  can  now  consider  multiphoton  pumping  upper 
vibrational  levels  of  simple  molecules,  exactly  on  resonance,  to  produce 
lasing  action  at  a variety  of  wavelengths  in  the  infrared. 

An  excellent  review  on  infrared  tunable  sources  and  current  spectro- 
scopy  research  is  given  by  Hinkley,  Nill  and  Blum.  A table  summarizing 
their  comparison  of  performance  characteristics  of  various  tunable  infrared 
sources  is  shown  in  Table  I. 

Important  materials  used  in  nonlinear  optics  devices  are  tabulated  in 
Table  II.  Of  the  approximately  forty  chalcopyrite  crystals  only  four  are 
used  for  nonlinear  optics:  ZnGeF£,  AgGaSg,  CdGeASg  -"nd  AgGaSe2-  The 
enormous  output  powers  and  high  efficiency  of  C0£  lasers,  coupled  with  the 
large  nonlinear  figure  of  merit  of  CdGeAs2,  make  this  particular  combination 
a very  promising  source  for  high  power,  tunable  infrared  radiation. 
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TABLE  I 


Source 


Semiconductor 
Diode  Laser 


j Spin-Flip 
Raman  Lasers 


Gas  Lasers 
Zeeman  Tuned 
High  Pressure 


Nonlinear 

Optics 


Wavelength 
Range  (^m) 


1-34 


3 - 9 

9 - 11 


0. 3 - 30  n 


Highest  _1  Output 

Resolution  (cm  ) Power 


100  W pulsed 


3 X 10" 
3 x 10" 


1 KW  pulsed 


0. 3 x 10 


1 - 10  mw  cw 
10-100  MW  pulse 


1 uw  - 1 MW 


TABLE  II 

NONLINEAR  OPTICAL  DEVICES 


Source 


Parametric 

Oscillator 


difference 

Frequency 

Generation 


Wavelength 
Range  (|im) 

Material 

Pump 

0.4  - 4 

Li 

Ruby 

0.4  - 3.6 

LiNbO^ 

Nd:Glass/YAG 

1.2  - 8.5 

Ag3AsS3 

Nd-.CaWO^ 

3.4  - 7-9 

CdSe 

Dy:CaF2 

9-8  - 10.4 

CdSe 

Nd:YAG 

3-0  - 4.5 

LiNb03 

Ruby/ dye 

3.2  - 5.6 

Ag.AsS 

Ruby/dye 

4.1  - 5-2 

LiI03 

Ruby/dye 

10.1  - 12.7 

Ag3AsS3 

Ruby/dye 

4.6  - 12.0 

*AgGaS2 

Ruby/ dye 

7.0  - 15.0 

*AgGaSe2 

Nd: YAG 

11.4  - 16.8 

*CdGeAs2 

co/co2 

Two  Photon 
Mixing 


9-0  - 11.0 


Four  Photon 
Mixing 


2.0  - 21.0 


K vapor 


Dye/dye 


Chalcopyrites 
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B.  CRYSTAL  GROWTH 


1.  CdGeAs 

During  this  program  we  have  evaluated  three  growth  techniques  for 
CdGeAs^ , vertical  Bridgeman  growth,  Czolchralski  growth,  and  growth  from 
Bismuth-Cd  solution.  Our  results  show  that  CdGeAs^  self  nucleates  on  the 
container  walls  and  grows  rapidly  in  a nearly  dendritic  way.  These  experi 

\ 

ments  suggest  that  serious  consideration  be  g^ven  to  an  oven  design  with  a 
very  steep  gradient  or  to  the  proper  design  of  a Czolchralski  growth  furnace 
for  CdGeAs^.  During  the  study  program  we  did  obtain  a single  crystal  ap- 
proximately 1 cm  X 1 cm  x 15  nun  in  size.  When  fabricated,  this  crystal 
yielded  an  oriented  crystal  6 mm  long  by  8 mm  x 8 mm  area.  We  are  now 
evaluating  this  crystal  for  CO^  laser  second  harmonic  generation  (SHG). 

The  enormous  nonlinearity  of  CdGeAs2  make  further  growth  studies  important 
if  future  advantage  is  to  be  taken  of  this  unique  material. 

2 . AgGaSe^  and  AgGaS,, 

In  a previous  progress  report  we  included  a preprint  of  the  paper 

"Growth  of  AgGaSe2  for  Infrared  Applications"  by  R.  K.  Route,  R.S.  Feigelson, 

29 

and  R.  J.  Raymakers . This  paper  described  the  growth  of  single  crystal 
AgGaSeg.  Since  then,  we  have  extended  our  work  to  include  the  growth  of 
AgGaS  and  have  successfully  grown  single  crystals  up  to  1 cm  diameter  by 
1 cm  in  length. 

During  our  previous  growth  work,  we  noted  the  consistent  poor  optical 
quality  of  both  AgGaSe0  and  AgGaS,,  due  to  optical  scattering  centers.  For 
example,  the  measured  optical  loss  at  1.15  pm  due  to  scattering  was  2 cm"1 
for  both  crystals. 

\ 
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In  an  attempt  tc  improve  the  optical  quality  we  tried  two  growth  modi- 
fications. The  first  was  to  dope  AgGaSe2  with  T£  in  order  to  improve 
the  growth  rate.  This  procedure  proved  not  to  be  of  value.  The  second 
was  to  rapidly  quench  the  crystals  from  a high  temperature  (~  JOO°C)  to 
room  temperature.  The  quenching  procedure  did  improve  the  optical  quality 
by  an  order  of  magnitude  reducing  the  loss  to  near  0.2  cm  — 0.1  cm  at 

1.15  pm. 

This  remarkable  improvement  in  crystal  quality  has  led  to  a careful 
study  of  the  quenching  process.  The  work  is  still  in  progress  but  the 
main  points  are  clear.  Both  AgGaSe^  and  AgGaS2  suffer  a phase  separation 
upon  slew  cooling  which  leads  to  optical  scatter  centers.  The  phase  sepa- 
ration occurs  near  700°C  for  AgGaSe^  Crystals  can  be  cycled  through  the 
quenching  process  and  reproduce  the  desired  optical  quality  improvement. 
And,  finally,  large  single  crystals  can  be  quenched  without  cracking. 

The  details  of  this  important  crystal  quality  breakthrough  are  being 
studied  and  prepared  for  publication. 


II.  THE  TEA-CO,,  LASER  SYSTEM 

A . INTRODUCTION 

The  literature  on  CC^  lasers,  both  conventional  low  pressure  and 
pulsed  atmospheric  TEA  discharges,  is  vast.  This  report  describes  the 
high  peak  power  short  pulse  laser  system  engineered  at  Stanford.  The 
description  of  the  1 nsec  system  is  more  from  a practical  construction 
point  of  view,  and  ample  reference  to  the  existing  literature  is  provided 
to  supplement  the  theoretical  aspects  of  the  first  system's  operation. 

A discussion  of  a 10  atmosphere  CO^  laser  designed  at  Stanford,  and  an 
alternative  scheme  to  laser  triggered  spark  gap  switching  is  given 
elsewhere 

B.  1 nsec  TEA-C0o  LASER -AMPLIFIER  SYSTEM 
1 . Background 

This  laser  system  was  engineered  and  constructed  with  two  objectives 
in  mind:  first,  to  study  higher  order  nonlinear  optical  processes  which 

become  very  efficient  when  high  peak  power  pulses  are  used.  However,  since 
all  materials  will  burn  at  approximately  1-4  j/cm  input  energy  density, 
it  is  necessary  to  reduce  the  laser  pulse  width  to  limit  the  total  energy. 
Secondly,  it  was  necessary  to  answer  the  many  questions  that  needed  investi- 
gation, particularly  in  producing  stable,  large  volume  gas  discharges  at 
high  pressures  and  in  producing  short  pulses  at  10.6  pm.  Laser  development 
in  this  field  in  the  past  four  years  has  been  impressive,  with  output  powers 
now  in  the  gigawatt  range  and  energies  of  500  J in  one  nanosecond  long  pulses. 


There  are  several  excellent  review  articles^1- ^3  that  describe  the 
development  of  C02  lasers  from  the  first  cw  laser  of  Patel^  to  the 
present  300  J system  in  Canada.^ 

The  basic  problem  of  providing  uniform  electrical  discharges  in  large 
volumes  of  high  pressure  (l  - 20  atmospheres)  gases  has  been  attacked  in 
many  ways.  Electrode  structures  for  transverse  excitation  range  from  pin 

discharge  geometries, ^ 38  three  electrode  structure s^"4^  an(j  resistive 
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electrodes,  to  more  exotic  structures,  including  electron  beam 
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controlled  discharges,  which  are  not  discussed  here.  An  alternative 
approach  to  electron  beam  pumping  for  producing  stable,  large  volume  dis- 
charges is  the  use  of  ultraviolet  radiation  to  preionize  the  laser  medium. 
Measurements  of  the  electrical  characteristics,^0  of  the  spatial  and  temporal 
behavior  of  the  gain,^  ^ and  of  the  spectral  content^  of  TEA-CCh,  dis- 
charges have  been  made  by  several  workers.  Finally,  several  authors  have 
attempted  to  fit  theoretical  models  to  TEA-C02  las.  r characteristics. 66“71 
Two  1 nsec  C02  laser  systems  of  comparable  design,  although  with  higher  output 
powers,  have  been  discussed  in  the  literature.72-^ 

2.  Design  and  Operation 

The  overall  system  design  is  shown  schematically  in  Fig.  3. 
a.  Laser  Oscillator 

The  pin  discharge  oscillator  is  similar  in  design  to  one  built  by 
Robinson.  It  consists  of  400  - 1,000  0 resistors  connected  in  parallel 
for  a cathode  and  an  aluminum  plate  for  an  anode.  We  obtained  a stable,  arc- 
free  oscillator  when  the  resistors  were  cut  off  with  equal  length  tips  and 
both  electrodes  are  mounted  on  a plexiglas  frame  to  keep  them  parallel.  The 
10-meter  curvature,  high  reflecting  mirror  is  internal  to  the  laser  tube. 


GENERAL  LAvOUT  OF  1 nsec  TEA  LASER  SYSTEM 
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FIGURE  3. 


General  layout  of  1 nsec  TEA  laser  system. 


A sodium  chloride  window  at  Brewster's  angle  provides  polarized  output 
and  an  adjustable  aperture  maintains  TEMqo  mode  operation.  The  output 
mirror  is  a flat  65$  reflecting  germanium  mirror  located  141.61  cm  from 
the  high  reflector.  Output  pulses  from  such  an  oscillator  are  typically 
200  nsec  in  duration  srd  0-5  - 1.0  MW  peak  power  in  TEM  q mode,  at  a 
repetition  rate  of  several  pulses  per  second. 

Since  the  oscillator  has  been  used  for  higher  order  nonlinear  optics 
processes  investigations,  we  wanted  as  stable  an  output  purse  as  possible. 

At  low  enough  gains,  namely  when  the  discharge  voltage  is  low  or  the  mirrors 
slightly  detuned,  no  mode  beating  occurs  and  this  oscillator  runs  an  approx- 
imately yja  peak- to- peak  stability  over  sixty  minutes  of  operation.  Figure  4 
shows  the  re  .suits  of  a measurement  of  the  laser  stability,  he  used  a 
crystal  of  CdGeAs^  to  produce  THG  of  the  laser  output.  The  third  harmonic 
signal  was  monitored  as  a function  of  time  and  shows  15$  peak-to-peak 
fluctuations.  Moreover,  by  doing  SHG  in  the  same  crystal  and  allowing  the 
fundamental  laser  frequency  to  mix  with  the  SHG  signal,  we  again  generate 
a 3 cu  signal.  Here,  however,  the  fluctuations  should  be  proportional  to 
twice  the  laser  fluctuations  rather  than  three  times,  as  in  the  THG  case. 

The  result  of  these  measurements  is  that  the  laser  operates  in  TEM  mode 

00 

at  5$  peak-to-peak  stability  over  long  periods  of  time.  At  higher  voltages, 
and  with  proper  alignment  mode  beating  occurs  and  the  laser  output  stability 
degrades  considerably.  Finally,  when  completely  modelocked,  the  pulse-to- 
pulse  stability  is  again  about  5$  . 

b.  Modelocking 

The  modelocker  consists  of  a germanium  acousto-optic  modulator^"’’ 
block  driven  by  a 36°  y-cut  LilfbO^  transducer  operating  as  a Bragg 


LASER  FLUCTUATIONS  =£>  5% 


deflector  = 41.4'  inside  and  2°  46'  outside  the  crystal).  Light 
enters  the  modulator  in  the  [110]  direction  polarized  along  the  [111] 
direction  for  maximum  interaction.  The  acoustic  driving  frequency 
50.313  MHz  is  half  the  cavity  mode  frequency  (e/2L)  , In  designing  such 
a modulator  with  high  index  of  refraction  material  (n  = 4.00  at  10.6  pm  in 
Ge),  one  must  consider  alternative  designs  to  Brewster  angle  cut  pieces 
(®B  = 76°  in  Ge)  in  order  to  make  reasonably  small  modulators.  By  making 
plane  parallel  faces  the  size  of  the  modulator  is  reduced  considerably  but 
etalon  effects  may  be  introduced  into  the  laser  cavity.  By  anti- reflect ion 
coating  the  modulator  faces,  the  bandwidth  of  this  etalon  can  be  made  several 
times  larger  than  that  of  the  laser  cavity77  and  therefore  reduce  the  problem 
of  cavity  instability.  The  same  modulator  with  no  anti- reflection  coating 
seriously  affects  laser  stability,  since  the  bandwidths  are  then  comparable. 
Figure  5 shows  the  pulse  train  detected  by  a fast  strontium  barium  niobate 
(SBN)  pyroelectric  detector78  and  displayed  on  a Tektronix  7904  oscilloscope. 

uPP®r  trace  shows  a superposition  of  10  modeloc  ced  pulse  trains, 
illustrating  both  very  clean  modelocking  and  small  peak-to-peak  fluctuations. 
The  lower  trace  shows  a 2 nsec  pulse  width  which  is  the  limiting  resolution 
of  the  7904  plug-in.  The  actual  pulses  are  1 nsec  wide. 

c.  Pulse  Selector 

The  train  of  1 nsec  pulses,  separated  by  10  nsec,  passes  through  a 
pilse  selector  consisting  of  a CdTe  electro-optic  crystal,79"82  a crossed 

QO  Ah 

polarizer,  laser  triggered  spark  gap,  and  coaxial  cable  charged  to 

the  correct  voltage.  The  modelocked  train  enters  the  electro-optic  crystal, 
r«^lects  off  the  Ge  Brewster  angle  plate  polarizer  and  triggers  the 
pressurized  spark  gap.  A length  of  coaxial  cable,  charged  to  twice  the 
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FIGURE  5 


(b) 

(a)  Modelocked  COg  pulse;  (b)  individual  pulses,  resolution 
limited  by  oscilloscope  preamp  bandwidth. 
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half  wave  voltage  of  the  electro-optic  crystal,  then  discharges  across  the 
crystal.  By  picking  the  length  and  impedance  of  the  cable  correctly  ('  i 
this  case  2 ft  of  73  fi  cable),  one  can  impress  a square  wave  of  approximately 
1 nsec  rise  and  fall  times  and  10  nsec  base  width  on  the  crystal.  The 
polarization  of  a single  1 nsec  pulse  is  then  changed  by  90°,  allowing  it 
to  pass  through  the  Ge  polarizer  and  out  to  the  C02  amplifier.  The 
CdTe  electro-optic  crystal  is  46  mm  long  and  5 X 5 mm  in  cross  section. 

The  half  wave  voltage  is  4.35  kV,  so  the  transmission  line  must  be  charged 
to  at  least  8.7  kV  (twice  what  appears  across  the  crystal).  Light  propagates 
in  the  [110]  direction  in  the  crystal,  with  both  electrical  pulse  and  light 
polarization  along  the  [110]  direction  for  maximum  interaction.®^  The  CdTe 
crystal  is  mounted  in  a brass  plate  electrode  structure  designed  to  match 
the  characteristic  73  n impedance  of  the  transmission  line.®5 

The  pulse  selector  assembly  is  by  far  the  most  unstable  part  of  this 
entire  system.  It  is  quite  difficult  to  trigger  the  spark  gap  reliably  over 
more  than  100  pulses  without  readjustment  of  the  laser  focus.  More  Importantly, 
however,  such  a fast  spark  discharge  radiates  very  high  frequency  noise 
(1  - 10  GHz)  which  makes  detection  of  the  switched  out  1 nsec  laser  pulse 
very  difficult  on  uny  fast  detector.  By  using  a slow  but  sensitive  pyro- 
electric detector  as  an  energy  monitor  we  have  observed  stable  switching  of 
pulses  from  the  modelccked  pulse  train.  By  replacing  the  spark  gap  with  a 
completely  electronic  detection  and  switching  system,  it  appears  possible  to 
make  a noise  free,  reliable  pulse  selector.  A suggested  design  is  discussed 
in  reference  30. 


d.  Amplifier 


The  single  1 nsec  pulse  then  passes  through  a three  pass,  modified 
Rogowski  profile  CC>2  amplifier.  To  make  the  discharge  more  uniform, 
two  arrays  of  spark  sources  called  "sliding  sparks"87"89  line  the  sides 
of  the  electrode  structure.  The  spark  sources  were  made  in  a transmission 
line  geometry  to  allow  easy  handling.  By  milling  a pattern  of  rectangular 
tabs  on  the  top  of  double  copper  clad  circuit  board,  one  creates  a trans- 
mission line  with  gaps,  which  can  be  modeled  as  an  array  of  capacitors 
coupled  together  by  the  st-ay  capacitance  of  each  gap.  A pulse  propagates 
on  this  line  by  arcinto  ^ross  each  gap,  thus  emitting  UV  radiation  into  the 
laser  plasma.  The  geometry  of  the  electrodes  for  both  laser  oscillator  and 
amplifier  is  shown  in  Fig.  6.  The  amplifier  has  two  internal  mirrors  made 
by  depositing  gold  coatings  on  flat  glass  so  that  the  amplifier  can  be  used 
for  1.  2,  3 or  5 passes.  We  have  measured  the  single  pass  gain  to  be  l#/cm 
of  excited  volume.  Topical  C02  amplifier  gains  range  from  l£/cm  to  1^/cm, 
depending  on  gas  mixture,  electrode  structure,  preionization  techniques,  and 
discharge  voltage.  In  this  amplifier  we  are  limited  in  voltage  to  30  kV  by 
our  power  supply  and  can  only  operate  at  300  torr  pressure.  The  measured 
gain  is  quite  good  under  these  circumstances.  It  would  be  necessary  to 

modify  the  high  voltage  electronics  to  allow  operation  at  higher  voltages 
and  pressures  to  raise  the  gain. 

In  contrast  to  the  resistively  ballasted  oscillator  which  operates  on 
C02  and  He  alone,  the  amplifier  requires  a 2 : 1 ; 3 ratio  of  C02,  and  He 
to  operate  stably.  Since  both  oscillator  and  amplifier  power  supplies  are 
the  same,  it  appears  that  the  impedances  of  each  electrode  structure  are  not 
the  same.  The  switching  thyratron  likes  to  look  into  a 50  n load;  apparently 
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PLEXIGLAS  TUBE 
AND  SUPPORTS 


AMPLIFIER  CROSS  SECTION 


FIGURE  6.  Oscillator  and  amplifier  electrode  structures 


by  adding  nitrogen  to  the  amplifier  medium  we  balance  the  capacitance  of  the 
electrode  structure  to  make  the  amplifier  look  more  like  a 50  ^ load. 


When  properly  impedance  matched,  both  devices  emit  relatively  little 
electrical  noise  by  comparison  to  the  laser  triggered  spark  gap. 

e.  Electronics 

The  circuit  diagrams  for  the  system  rate  generator  and  master  control, 
thyratron  trigger  circuits,  and  high  voltage  electronics  are  shown  in 
Figs.  7 to  9-  The  modelocking  electronics  is  not  shown  here.  It  consists 
of  a standard  Tektronix  190B  signal  generator  operating  at  25-157  MHz,  an(i  a 
Heathkit  DX-60  transmitter  modified  to  operate  in  pulse  mode  at  twice  the 
oscillator  frequency.  Since  the  effective  impedance  of  the  acousto-optic 
modelocker  is  2.2  fi,  it  is  necessary  to  use  a pi  network  to  match  the 
50  n output  impedance  of  the  Heathkit  transmitter. 

Since  we  require  3.7  kV  on  the  transmission  line  for  operation  of  the 
CdTe  electro-optic  crystal,  and  we  must  maintain  a rise  time  of  1-2  nsec, 
it  seemed  impossible  at  first  to  use  an  avalanche  transistor  array  for  the 
pulse  selector  electronics.  Consequently  we  had  to  use  a laser  triggered 
spark  gap  as  our  fast  switch.  DC  voltage  is  provided  by  a Variac  and  DC 
power  pack.  However,  by  using  a transmission  line  transformer  it  now  appears 
that  an  array  of  approximately  32  avalanche  transistors  can  be  used  to  generate 
4-5  kV  across  the  electro-optic  crystal.  A block  diagram  of  the  proposed  1 

electronics  is  shown  in  reference  30* 
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nsec  TEA  LASER 


FIGURE  7.  Master  control  circuit  for  1 nsec  laser  system 
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HIGH  VOLTAGE  ELECTRONIC 


III.  NONLINEAR  OPTICS  THEORY 


This  section  presents  a theoretical  treatment  of  higher  order  non- 
linear optical  processes  in  crystals  witn  h2  m symmetry.  It  begins  by 
determining  the  equations  for  second  ( SHG ) , third  (THG)  and  fourth  (FHG) 
harmonic  generation  in  a crystal  with  absorption  at  the  fundamental  fre- 
quency co  and  each  of  its  harmonics  2cd,  3cd  and  )+cd  . Also  considered 
is  the  case  of  nonphasematched  mixing  to  produce  higher  harmonics,  namely 

CO  + 2m  - 3cd  , CD  + 3co  ->  IkD  , and  2cd  + 2cd  - ku  . All  equations  are 
derived  in  cgs  units. 

A.  SECOND,  THIRD,  AND  FOURTH  HARMONIC  PROCESSES 

1 • Nonlinear  Polarization  Equations 

Throughout,  this  discussion  follows  the  works  of  Midwinter  and  Warner 
on  second  order?0  and  third  order91  ptoccsses>  keepi„g  ^ ^ ^ ^ 

are  missing  a factor  of  2 for  the  case  of  second  order  mixing,  and  factors 
of  or  6 for  third  order  mixing. ^ 

We  can  write  the  general  tensor  equations  for  this  crystal  class  for 
second,  third  and  fourth  harmonic  generation  processes  as 


P,  = d.  E.  E, 

ijk  j k 


Pi  Cijk£  Ej  Ek  Ei 


(3-la) 

(3-lb) 


>.  ..  „ E.  E,  E E 
J k l r 


(3.1c) 


If  we  restrict  ourselves  to  frequencies  in  the  transparency  range  of 
the  material  and  assume  Kleinman  symmetry,^  the  nonzero  elements  of 
the  d,  c,  and  b tensors  can  be  determined  and  are  shown  in  Fig.  10. 

The  reduced  subscript  notation  is  explained  in  Fig.  11  for  ell  three 
processes.  In  Appendix  I a practical  discussion  is  presented  for  cal- 
culating higher  order  nonlinear  coefficient  tensors  from  those  of  lower 
rank.  Equations  (3.1)  are  written  out  in  detail  as  shown  in  Fig.  12. 

The  cases  of  SHG  and  THG  have  been  discussed  by  Kildal^  and  are  simply 
summarized  here.  The  case  of  FHG  is  treated  in  detail. 

2.  Fourth  Harmonic  Generation 
a.  Effective  Coefficients 

With  the  exception  of  Se,  CdGeAs2  is  the  only  crystal  that  phase- 
matches  for  FHG  in  the  infrared.  From  Fig.  10  we  see  that  there  are  two 
independent,  nonzero  components  of  the  FHG  tensor  for  a crystal  with  42  m 
symmetry,  assuming  Kleinman  symmetry.  For  convenience,  let  us  define 

bl  (*0  = ba  bl(5)*bB 

Then,  from  Figs.  10  and  12(b)  we  have 

Px  = ba  tVz1  + b9  !UEzEy  + 12  VVv 

Py  - *>„  + bB  [ks^l  * 12  ExBzEp  (3. 3) 

Pz  * ba  [15EzEy^!  * VbExE)  + ltEyE^ 
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FIGURE  10.  Nonlinear  coefficient  tensors  for  2u>,  ya  and  lw>  processes  in  k2m  crystals. 


FIGURE  11.  Reduced  subscript  notation  for  2cjj,  iko  tensors 


For  FHG  we  have  four  possible  types  of  phasematching  situations.  They 
are  given  below  with  the  corresponding  index  matching  equations  for  a 
positive  birefringent  crystal. 


T^ype  I e + e + e+  e-»o 
Type  II  e+e  + e + o-*o 
Type  III  e + e + o + c->o 
Type  IV  e + o + o + 0-»o 


4n  (6,0))  = 4n  (2ko) 
c o 

3ne(fi,o>)  + nQ(o))  = 4n  (iku) 

° (3-4) 

2ne(A,a)  + 2nQ(m)  = 4n  (1«d) 
ne(9*u>)  + 3n0(a>)  = 4nQ(4a)) 


Here  e and  o refer  to  extraordinary  and  ordinary  waves.  Only  the 
first  two  types  phasematch  in  CdGcASg. 

We  must  now  calculate  the  effective  nonlinear  coefficients  for  the 

Type  I and  Type  II  processes.  From  Eq.  0.4)  for  the  Type  I process  it 
is  evident  that 


P(4<u)  = Po(lku)  = 


P(4<c) 


E(o)  = Ee(oo)  = 


E(o>) 


(3-5) 


The  geometry  is  explained  in  Fig.  1?.  Combining  Eqs.  (3.3)  and  (3.5)  we  have 


P>: 


F - - 4basin-e  cos  fisin  0 - 4b3[cos36sin  Qsiri3£75  + 3 cos39sin  6cos20  sin  0] 


,3„_* 


E 

P 

E 


""1+  ~ 4b^sin  6 cos  0cos  0 - lb g[cos3Bsin  0cos30  + 


3 2 

3 cos  fisin  Psin  0 cos  0] 


(5.6) 
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we  have 


P = 


P sin  0 - 
x r 


V cos  0 

y 


P sin  0 


= - [4b  sin^ecos  ft  + 4b„cos^ftsin  0(sin^  0 + 3 0)]  sin^  0 


. 2 


E 


dv*vu  w • -r*/^ 


P cos  0 - O 2 2vvn2w( 

JL_ = [4b  sin^ftcos  6 + 4bQcos-39sin  ft(3  sin  0 + cos  0)  J cos  0 

E4  “ P 


(3*7) 


We 


define  the  FHG  effective  coefficient  as  P0(4cn)  = beff.  E (o>)  and  from 


Eq.  (3-7)  we  determine  it  to  be 


beff(I)=  S[b 


a 


■ 2 . 
sin  0 


O 

+ bQ  cos  0]  sin  2 0 cos  2 0 

p 


(3-8) 


for  Type  I phasematching  in  crystals  with  42  m symmetry.  It  should  be 
noted  that  for  43  m symmetry  this  expression  reduces  to  bef£.(43  m)  = 
2ba  sin  2 6 cos  2 0.  In  a similar  way  we  calculate  the  Type  II  effective 
coefficient  to  be 


n ,,  „ 

b (ll )=  4[b  sin  9 + b cosc  ft]  sin  0 sin  2 0[sin^  pcos  p] 


eff 


(3-9) 


The  angles  0,  0,  and  p are  again  defined  by  Fig.  13  and  E = E sin  p' 


,-JS 


It  is  easy  to  show  that  the  expression  sin^  p cos  p is  maximized  when 


„f~ 

p = + + nn  , in  which  case  4 sin^  p cos  p = . Moreover,  cos  2 0 

is  maximum  for  0 = 0,  whereas  sin  2 0 is  maximum  for  0 = r 


Therefore 


the  maximized  equations  in  terms  of  9 alone  are 


2 4 

Type  I P(4oi)  = 2[ba  + (bp  - bQ)  cos  9]  sin  2 6 E (co) 


Type  II  P(4<u)  = ^ [ba  + (bp  - bQ)  cos2  P]  sin  9 M 


(3.10) 


From  index  of  refraction  data  due  to  Boyd'*"'1  we  calculate  a phasematching 
angle  for  \ = 10.6  pm  of  ®m  = 67.85°  in  which  case 

Type  I P(4<u)  = 1.406  [bQ  + .l42(b^  - ba)]  E^(oo) 


Type  I P(4co)  « 1.406  ba  E (a>) 


(3*11) 


for  x = 10.6  pm  . Two  problems  arise  in  looking  at  this  phasematching 
process.  There  is  a good  possibility  of  wave  vector  walk-off  in  any  crystal 
long  enough  to  give  appreciable  signal.  One  can  pdssibly  generate  4o>  by 

* 

mixing  3co  + to  -»  tao  since  THG  is  phasematched  theoretically  at  64°  for 
X = 10.6  pm  and  the  two  processes  could  overlap.  From  Boyd's  index  data 
the  turning  point  of  the  Type  II  FHG  curve  occurs  at  X = 13  pm  and  P^  = 82°. 
In  this  case 


Type  II  P(4cu)  « 1.29  b„  E (oj) 


(3.12) 


for  x = 13  n®.  Even  though  the  Type  II  process  has  the  advantage  of  being 
phasematched  close  to  90°  so  that  walk-off  is  not  a problem,  it  is  nevertheless 
quite  difficult  at  present  to  consider  generating  enough  radiation  at  13  pm 
to  actually  observe  the  Type  II  - FHG  process. 
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Since  the  index  of  refraction  measurements  for  this  crystal  are  not 


yet  that  accurate,  we  must  consider  these  FHG  calculations  as  approximate. 
However,  the  fortunate  proximity  of  phasematched  Type  II  - THG  with 
Type  I - FHG  bears  further  investigation.  In  the  case  of  producing  ya 

7 

by  mixing  2cjo  + 03  3^  > the  THG  output  depended  on  the  oHG  signal. 

Similarly,  here  it  is  possible  to  consider  the  mixing  process  30  + cd  -»  *kD 
where  the  FHG  output  will  depend  on  the  THG  signal.  We  have  in  fact  observed 
FHG  due  to  this  mixing  process. 

b.  FHG  Fower  and  Conversion  Efficiency 

We  can  calculate  the  FHG  power  at  phasematching,  taking  into  account 
the  properties  of  a focused  beam  and  the  transmission  of  the  crystal  faces. 
Since  for  all  crystal  sizes  of  interest  here  the  effective  length  of  the 
focal  region  is  much  larger  than  the  crystal  length,  we  can  use  a plane  wave 
analysis.  The  electric  field  in  a medium  can  be  written  in  this 

approximation  as 

E(oj)  = A(z,cjq)  exp (~A  Q^z)  expti^z  - cot)]  (3-13) 

for  propagation  along  some  direction  z.  The  FHG  polarization  is  given  then 
by 

P^Clkn)  = beff  A^(gd)  exp(-2a1Z)  exptiClk^z  - Uot)]  (3.14) 

In  the  thin  slab  between  z'  and  z'+  d z',  the  increment  of  the  FHG  signal 

9U-97 

can  be  obtained  from  Maxwell's  equations  and  is  given  by 


2iti 


dA(z,4oj)  = 


r%> 


exp(-ik,  z'  + i 4<d  t)dz'  (3*15) 


where  the  geometry  is  shown  in  Fig.  I1*  and  n^  is  the  refractive  index 
at  = Iko  . This  generated  signal  is  somewhat  absorbed  in  the  (£  - z'  ) 
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distance  remaining  in  the  material.  Consequently  the  increment  in  the 
FHG  signal  at  the  end  of  the  crystal,  due  to  that  produced  in  the  thin 
slab,  is  given  by 


^jti  a. 

dA(z,4a>)  = r (4to)  exp(-ik,  z'  )exp[-|a,  (£-z'  ) ]dz'  exp(ikdt)  (3. 16) 


\ C 


2iti  1124 

dA(z,Iko)  = b ff  A (to)  exp(-ia^)  exp(-Aaz' ) exp(iAkz' )dz'exp(i4a>t) 

. n.  c 


where 


(3.17) 


Aa  = 2a±  - and  Ak  = 4^  - k^ 


Ignoring  the  time  dependence  we  have 


4.  2jti  a>|^  , £ 

A(/,4<d)  = J dA(lku)  = beff  A (cD)exp(-|a^f ) f exp(iAk-Aa)z'dz 

0 L n.  C J o 


(3-18) 


Integrating,  we  get 

[ 2ni  cdv  . 

A(£,kiu)  = I b . A (oj)  exp(-| 


n4  c 


exp(iAk  - Aa).g  - 1 
(iAk  - Aa) 


(3.19) 


|E(£,4o))|2  = j A(i,4cu)| 2 = A(£,hw)A*(£,bw) 


, v . 2 J2-4  4,  l2  fcosh(Aae)  - cos(Aki) 

|E(/,4<d)|  = 2 beff  A (ai)  exp  (-Of.  i)exp(-AoX)  

V (Ak  f + (Aa)2 


We  can  write  the  power  generated  at  Iko  as  P(4m)=A  < S > where  S is 

JtW? 

the  Boynting  vector  S = ^ j E| 2 and  A = — ^ for  a Gaussian  beam. 

2 


PCM  = \ lECto)^  = ^ *\  |a(1m>)| 


(3.21) 


where  w^  is  the  beam  waist.  Since  the  electric  field  amplitudes  are 
related  by  transmission  factors  at  the  faces  we  can  write 


(3-22) 

as  the  general  expression  for  the  FHG  power  coming  out  of  a crystal.  We 

have  used  the  facts  that  E(o>)  (inside)  = E(od)  (outside)  and 

E(1ks)  (outside)  = — e(Uo)  (inside)98  and  that  - (Aa  + ak ) = - (2a  + A a ). 

If  we  consider  the  case  of  phasematching,  Ak  0 , then  we  can  simplify  the 
last  term  as  follows 


cosh(Aki)  - cos(Aki) 
(Ak)2  + (toe  f 


1 + A(Aai)2  - l + |(Ak if 


(Ak)2  + (Aa)2 


« l2 


(3-23) 


I 


Therefore  for  phaseraatched  FHG  we  have 


(3-25) 

where  ^ + n and  ^ + n from  the  Fresnel  equations,58  and 

n4  = n for  phasematching.  The  FHG  conversion  efficiency  is  given  by 


(3-26) 


c.  Angular  Half  Width 

4 

Kildal  has  given  a clear  description  of  the  calculation  for  the 
angular  half  width  of  generated  SHG  and  THG  signals.  Such  a measurement 
is  important  for  several  reasons.  First  it  serves  as  a check  on  the  type 
of  process  one  is  actually  observing;  mixing  and  harmonic  generation,  for 
example,  have  different  angular  tolerances.  Second,  it  is  an  indication 
of  the  practical  angular  tolerance  limits  one  encounters  in  turning  a 
crystal  to  phasematch  different  input  frequencies.  Third,  it  serves  as  one 


I 


98 

and  the  Foynting  vector  walk- off  angle  T}  as 


tan  T|  = 


n*(e)  /I 


n 


sin  2 6 


n 


(3.30) 


From  these  two  expressions  one  finds 


bng(0 ) 

= - n (6)  tan  T] 

be 


(3.31) 


Finally,  for  near  normal  incidence  A0  external  = nA0  internal.  Combining 
this  result  with  Eqs.  (3.28)  and  (3*31)  we  find 


X 

A0T  ext  " rs 

1 ext  8.9&g  tan  T| 

X 

n (M 

A0II  ext  “ c . 

6.73^  ^an  T] 

ne(cu,0 ) 

(3.32) 


3.  Summary  of  Harmonic  Generation  Results 

L L 

At  this  point  we  can  tabulate  the  results  for  SHG,  THG,  and  FHJ 
in  a crystal  with  \2  m symmetry.  The  effective  nonlinear  coefficients 
for  the  various  phf sematching  conditions  are  given  in  Table  3.  The 
expressions  for  output  powers  at  the  harmonic  frequencies,  in  terms  of 
the  input  power  at  00  , are  shown  in  Table  4.  Finally  the  angular  half 
widths  are  given  in  Table  5.  It  should  be  noted  that  the  widths  get 
narrower  as  we  go  to  higher  order  harmonics.  The  tuning  curves  for  these 
three  processes  are  shown  in  Figs.  15  and  16. 
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SHG  PHASE  MATCHING  ANGLE  (deg) 

FIGURE  15.  SHG  angle  tuning  curves  in  CdGeAs  . 
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IGURE  1 6.  THG  and  FHG  angle  tuning  curves  in  CdGeAs, 


SHG 


I 


d1  = cL^  sin2ecos  2 0 
II  l^ii  = sin  9 sin  20  sin  2p 


THG 

I 

cT  = - £ (cu  - 3 cl8)  cos3  e sin  4 0 

II 

CII  = [l  (cll  ‘ 3 C18}  C°s2  9 Sin2  2 ^ + Cl6  Sin2  9 

2 2 
+ c q cos  6 ] cos  p sin  p 

III 

III  = ¥ (cn  ‘ 3 cl8)  cos  6 sin  1+0  sin  p c°s2  p 

FHG 

I 

bj  = 2[b^(j+ ) sin2  9 + cos2  9]  sin  2 P cos  2 0 

II 

O Q Q 

bII  = ^^bi(4)sin  9+  bl(5)  cos  sin  9 sin2^  tsin  Pcos  p^ 

; 


TABLE  IV 


POWER  EQUATIONS  FOR  HARMONICS 


/l6VWp\/  2\/2tt  d i\  / \2 

p(a°)  -(-H  16  - )\ — — ) 

\c  AWj^/X  / \ n \ / \x  / 


t2)  P U)  exp[-  |(a2  + 2ax)  j] 


P^to)  exp[-  -|(a  + 30^}  i] 


P (cu)  exp[-  ^(a,  + ha.)  i] 


where  wx  -»  w^  are  beam  waists,  t±  -*  are  electric  field  transmission 
factors  and  -*  0^  are  loss  constants  for  beams  with  a>  _»  ko>  . & is 

the  effective  interaction  length. 


wauw j iy wjwwi.mui  .i**«ip,i.i.j mm.rp . *>■*■■  ■»-■* 

< 

f 

B.  MIXING  PROCESSES 

1.  a + 2cd  ->  3a> 

There  are  two  mixing  processes  that  can  generate  yn  and  4<o  which 
should  be  discussed.  The  first  process 

ca  + 2oj  -*  ya  (3*33) 


is  not  phasematched  and  occurs  along  the  direction  of  phasematched  SHG. 

By  using  two  sources  one  at  10.6  (im  and  the  other  at  5.3  (m  we  can  have 
phasematched  mixing.  In  terms  of  ordinary  and  extraordinary  polarizations 
the  only  combination  of  importance  here  it 

•1  U — \ 

Ee(ci3)  + Eq(2cjd)  -»  Eo(3cd)  • f3.3^) 

Because  of  positive  birefringence,  the  harmonics  must  be  ordinary  waves. 
Moreover,  ov.U.y  the  extraordinary  wave  at  cu  will  generate  a mixing  signal, 
since  Eo(co)  + Eo(2cu)  cannot  give  Eo(3oi).  Since  this  is  a second  order 

process,  we  have  some  effective  mixing  coefficient  dmix  (effective).  We 
can  determine  its  form  by  writing  down  the  expressions  in  Figs.  10-12 
for  a mixing  process. 


Px(3cd)  = 2dl4  [Ey(oj)  Ez(2o))  + Ey(2o>)  Ez(®)] 

Py(3ca)  = 2du  rEx(®)  Ez(2cd)  + Ex(2cd)  E^cu)] 

Pz(3<u)  = 2d11+  [Ex(«j)  Ey(2cD)  + Ex(2cu)  Ey(®)]  (3*35) 


Note  here  P(3o>)  = P0(3oo)  so  P (3®)  = 0 . Also  E(2cu) 
Ez(2aj)  = 0 . These  equations  reduce  then  to 

* 2dl%  Ee  >’*>  Eo,y<2“> 

'./I  - 2ai4  Ee,z<»’e>  V(2mi 
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E0(2a')  so  that 


(3-36) 


where  P (3o>)  = sin  0 p v(3a3)-cos  0 P0  (3cu)and  E (oi,6 ) = sin  p sin  6 E(co). 
ow  ' r o,x  r ° >y 

96 

Note  the  degeneracy  factor  2 is  in  agreement  with  other  authors.  Now 
from  Table  3 we  have 

p (2oj)  = - d^  sin  6 sin  2 0 sin  2 p E (ca)  (3-37) 

96 

and  we  can  calculate  the  magnitude  of  A(2cd)  to  be 


A(z'  2cd)  = i dSHG  A^  exp[-  |(ag  + 20^)  z' ] z'  (3-38) 


n2c 


where  dgHG  = sin  0 sin  2 0 sin  2 p. 

mixing  process  we  have  a coherence  length 


Now  we  assume  that  for  the 

i , where 
c * 


6 it 

Ak  = k + k - k = — 
3 X 


ni  + 2n2 


n„ 


X 

It 

c 

Ak 

6 

+ 2n2 

l 3 ' 

(3-39) 


and  that  the  magnitude  of  E(2cd)  is  constant  and  given  by  its  value  at 
(£  - | £c),  namely 

. 2 it  , 

A(£  - | £ , 2co)  = i d A (o>)  [£  - ± exp  [-  -(«2  + 20^)  £] 

n2c 

(3-40) 
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Following  the  same  procedure  as  for  the  FHG  case,  we  calculate  the  power 


at  (30)  due  to  mixing  to  be 


P(3o) 


/,*v2  , 2 > , 

/ 16  \ / w^  \ / 
\ c / \w?/  \ 


^IX  dSHG 
2 

n2n^ 


fc(f-2ic)  (tit3) 


exp  [-  |(«2  + 30^)  i] 


where  dgHg  ■=  - d^4  sin  q sin  2 0 sin  2 p 
dMIX  = " 2dl4  Sin  9sin  p 


2.  to  + 3ta  -»  4<p 

The  case  of  non-phasematched  mixing  of  a THG  signal  with  the  funda- 
mental to  produce  a signal  at  can  be  analyzed  similarly.  By  analogy 

to  Eq.  (3.34)  we  have 


Ee(cn)  + Eo(3o>)  -*  Eo(4o) 


The  coherence  length  for  this  process  is  given  by 


£c  ~ ->  -» 

* kx  + k3  - k4 


n , + 3 n - 
e,l  J 0.3 


Here  we  are  concerned  with  mixing  against  a signal  at  3a>  generated  in 
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a type  II-THG  process. 


The  calculated  output  power  at  4u>  is  given  by 


CTHG  = l[(cil"3Cl8)sin2  ^ C°s2fl+  Cl6  Sin2f5+  °l8  C°s2e]  Sin2  P C0S  P 
At  this  point  it  should  be  noted  that  the  process  2ai  + 2a3  ->  4oi  does  not 

work  because  all  three  waves  are  ordinary  waves  when  we  consider  mixing  of  a 
SHG  signal  with  itself  to  give  FHG.  On  the  other  hand,  if  we  have  an  external 

source  for  the  2<u  signals  we  can  polarize  the  input  waves  to  make  this 

mixing  process  possible  as  well. 

3.  Angular  Half  Widths 

We  calculate  the  angular  half  widths  for  these  two  mixing  processes 
by  the  same  procedure  used  in  Eqs.  (3.28)  to  (3*32) 


cn  + 2oi  -♦  3cs 


and 


0) 


nQ(3o3) 

11  ext  2.25  l tan  p 

ne(fi,a>) 

4o) 

(3.45) 


(3>6) 


wmm 
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IV.  EXPERIMENTAL  AND  THEORETICAL  WORK  ON  HIGHER  ORDER  NONLINEAR 

PROCESSES  IN  CdGeAs2 

A . INTRODUCTION 

Since  CdGeAs^  has  the  largest  figure  of  merit  for  second  order 
nonlinear  optics  processes,**  and  since,  as  shown  in  Section  III,  it  can 
be  pha sematched  for  third  and  fourth  harmonic  generation  processes,  we 
undertook  an  experimental  and  theoretical  investigation  of  these  two 
interactions.  The  results  on  THG  experiments  in  CdGeAs2,  and  the  theoretical 
models  used  to  calculate  contributions  to  the  third  order  susceptibility 
from  bound  electrons  and  free  carriers  are  discussed  in  Appendix  II. 

This  section  gives  a summary  of  the  application  of  the  Bond  Charge  Model^-101 
to  the  calculation  of  fourth  harmonic  tensor  coefficients.  The  results  of 
the  third  harmonic  investigation  of  CdGeAs2  are  discussed  and  the  potential 
power  outputs  at  the  third  and  fourth  harmonics  of  10.6  pm  are  examined. 


■ B.  THIRD  HARMONIC  EXPERIMENTS  AND  THEORY 

1.  Separation  of  Components  by  Symmetry  Arguments 

One  of  the  principal  points  made  in  Appendix  II  is  that  nonlinear 
aisceptibilities  of  various  orders  arise  from  different  mechanisms  in  a 
semiconductor  material.  In  particular,  we  note  that  isotropic  materials 
1 have  non- vanishing  odd  order  susceptibility  tensors  while  even  order  sus- 

ceptibilities equal  zero  from  symmetry  considerations.  In  particular,  for 
semiconductor  materials  odd  order  susceptibilities  arise  from  two  sources: 

! electrons  localized  in  crystal  bonds,  which  reflect  the  symmetry  properties 
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of  the  crystal  lattice,  and  the  free  carrier  cloud,  which  has  isotropic 
symmetry.  For  crystals  with  42  m symmetry,  like  CdGeAs^,  we  can  compare 
nonzero  third  order  tensor  elements  with  those  of  systems  with  other 
symmetries.  Such  a comparison  is  shown  in  Fig.  17  where  the  nonzero  third 
order  tensor  elements  for  uniaxial  h2  m symmetry  are  obtained  from  Fig.  10. 

The  fact  that  the  c^g  element  is  common  to  the  case  of  42  m and  isotropic 
symmetry  suggest  that  this  tensor  component  arises  from  a component  of  the 
crystal  with  isotropic  symmetry,  such  as  the  free  carrier  cloud.  The 
element  c^  - 3c^g>  however,  appears  to  arise  solely  from  a component 
reflecting  the  symmetry  properties  of  the  particular  crystal  class,  i.e. 
the  electrons  localized  in  the  bonds.  From  Table  3,  moreover,  it  is 
apparent  that  one  can  separate  these  tensor  components,  and  hence  the  bond 
electron  and  free  carrier  effects,  by  doing  different  phasematched  THG 
experiments.  These  points  are  summarized  in  Fig.  18.  As  discussed  in 
Appendix  II,  we  have  performed  phasematched  Type  I ard  Type  II  THG  to 
measure  the  third  order  tensor  coefficients.  Until  now  the  only  measurements 
of  third  order  susceptibilities  in  semiconductors  have  been  carried  out  by 
mixing  experiments  (2u^  -cjo^  = cjd^)  in  a nonphasematched  configuration. 
Consequently  it  was  not  possible  to  measure  the  tensor  coefficients  independ- 
ently, as  we  have  done  in  CdGeAs^  . 

2.  Free  Carrier  Effects  for  Third  Harmonic  Generation 

In  Section  III  of  Appendix  II  contributions  to  the  third  order  suscept- 

lO^-ioo 

ibility  from  free  carriers  are  discussed.  1 ' Since  we  were  doing  harmonic 

generation  rather  than  mixing  experiments  (where  frequency  differences  can 
get  close  to  inverse  relaxation  times),  the  only  contribution  to  the  nonlinear 
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FIGURE  17.  Analysis  of  THG  tensor  components  by  symmetry  arguments. 
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SEPARATION  OF  BOND  ELECTRON  AND  FREE  CARRIER  EFFECTS 

BY  PHASE  MATCHED  THG 

cim  =*  c?m  (BOND  ELECTRONS)  + cj^  (FREE  CARRIERS) 

cfmE  -=»  SYMMETRY  OF  CRYSTAL  CLASS  42m 

cFC  ^ SYMMETRY  OF  I SOTROPIC  ELECTRON  CLOUD 

TYPE  I THG: 

ceff  * 

TYPE  II  THG: 

ceff  * 

+ <Ci6>BE  + FC 


( c 1 1 ■ 3c i q ) 


18 


BE 


,cll  ' 3c18) 


BE 


+ (c16  ' c18) 


BE  + FC 


FIGURE  18.  Bond  electron  and  free  carrier  contributions  to  third  order 
susceptibilities . 
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susceptibility  comes  from  the  nonparabolic ity  of  the  energy  bands. 


We  have  extended  the  work  of  Kane110  and  Kildal4  to  the  calculation 


of  third  order  nonlinear  susceptibilities  in  ternary  chalcopyrites.  The 
Kane  model  has  been  used  successfully  to  describe  the  third  order  properties 
of  elementary  and  binary  compounds.  108'1;L5  Details  of  the  calculation  of 
the  energy  band  shapes  by  our  extension  of  Kane's  model  are  given  in  Appendix 


II.  It  has  been  demonstrated  by  several  authors  that  the  third  order 


susceptibility  due  to  nonparabolic  energy  bands  is  given  by107"10^ 


Cijki(uWV  = ± 


bSl 


24^70^0^(0^+  o^+u^)  k 


c,v 


E f (E  ) 

° c’v  r>k.bk.i>k  bk, 


i J k £ 


(4.1) 


where  fQ(E)  if  the  Fermi-Dirac  distribution  function.  Here  we  pick  the 
+ Sign  and  Ec  for  material,  and  the  - sign  and  Ey  for  p-type 

material.  A comparison  of  the  energy  band  diagrams  for  chalcopyrite.  and 
sphalerite  materials  is  shown  in  Fig.  19.  The  parameters  used  in  this  model 
for  CdGeAs2  are  also  listed.  As  discussed  in  Appendix  II,  the  results  of 
this  free  carrier  model  for  p-type  CdGeAs2  are 


C16^  = -2.5xio"27n 


- 4P)  - 2 X ^7  N 


C16}  = L5X10'27  N 


C18)  ' C16^  = 6'3  x lo”28  N 


(4.2) 


where  N is  the  free  carrier  density.  We  use  the  fact  that  for 

an  isotropic  medium,  in  this  case  the  free  carrier  cloud,  c = Jo  . 

11  18 


We  derive  the  following  values  for  c and  r *»,  -p 

IV  16*  ana-  rrom  the  free  carnei 


— j-w  — 

acdel  using  a carrier  concentration  of  li  . 1.21,  x lo16  holes/cs,3.  (note: 
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SPHALERITE  COMPOUND  CHALCOPYRITE  COMPOUND 


FOR  CdGeAs2  E6=0.53eV,  A=0.38eV,  S = -0.25eV 

E,  =0.20eV,  Ej=  -0.32eV 


FIGURE  19.  Band  structure  of  sphalerite  and  chalcopyrite  compounds. 


o 

the  value  N = 5 x 10  holes/car 
c1 ^ = -1.4-9  X 10  ^ esu 
c^g  = -0.5  X 10  U esu 
cl6  = -3.1  X 10" 11  esu 


given  in  Appendix  II  is  incorrect). 

Free  carriers 
(N  = 1.24  X 1016 

O 

holes/cnr ) 


These  should  be  compared  with  the  results  of  the  Bond  Charge  Model 
c = 1.8  x 10  ^ esu 

C3j8  = 3 X 10  ^ esu  Bond  Electrons 

cl6  = 1.4  x 10" 11  esu 

3-  Comparison  of  Theory  and  Experiment  for  Third  Harmonic  Generation 
We  can  compare  the  values  for  the  effective  c's  for  Type  I and 
Type  II  THG  by  using  the  expressions  in  Table  3,  keeping  in  mind  that  we 
must  use  the  free  carrier  and  bond  electron  values  in  accordance  with 
Fig.  20. 

^ 1 : Cj  = - £ (cn  - 3cl8)BE  cos39i 

Type  II  : cI];  = sin2  p cos  p coa2*^ 

+(ci6  ci8  cos2en)BE  * FC] 

(M) 

In  our  experiments  the  phasematching  angles  are  = 50°,  0 = 60°,  and 

cos  p = . 

ST 


- 58  - 


Cm 


The  crystals  used  were  cut  for  maximized  0 angle  in  each  case. 

Therefore, 

c^.  = -0.14  X 10  '*‘1  esu 

Calculated 

CIX  = -0.53  X 10" 11  esu 

(4.4) 

By  using  the  experimental  setup  shown  in  Fig.  20,  we  observe  both  Type  I 
and  Type  II  THG  signals  in  CdGeAs^.  As  described  in  Appendix  II  we  observe 
a signal  at  3to  by  nonphasematched  mixing  of  the  second  harmonic  signal  with 
the  fundamental  2oo  + co  = 3^  • Since  the  value  for  di4  of  CdGeAs2  is 
known  to  be  8.4  X 10  ^ esu,  by  measuring  direct  Type  II  THG  in  the  same 
crystal,  we  obtain  a value  for  the  ratio  c /d^.  Finally,  by  observing 
Type  II  THG  in  one  crystal  and  Type  I THG  in  a second  crystal  with  the 
same  experimental  setup,  one  can  get  the  ratio  cj/cjj-  A phasematching 
curve  for  Type  II  THG  and  its  P^(o>)  power  dependence  are  shown  in  Fig.  21. 
From  these  experiments  we  have  determined  and  cTI  to  be 

| c jj  < 0.6  X 10  esu 

Measured 

|c  I = 12.5  X 10  ^ esu 

11  (^-5) 


! 


Since  our  signal  to  noise  ratio  was  only  2:1  in  the  c^  measurement 
we  must  put  an  error  range  of  + 50$  on  this  measurement.  It  does,  nevertheless, 
compare  favorably  with  the  value  given  in  Eq.  (4.4)  from  the  Bend  Charge  Model 
calculation.  However,  our  measured  and  calculated  c^^  differ  by  more  than 
one  order  of  magnitude.  There  are  several  reasons  to  explain  this  discrepancy, 
all  of  which  bear  further  investigation. 
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TYPE  II  THG  TUNING  CURVE 


■•(5) 


Pr  owl 

FIGURE  21.  THG  tuning  curve  and  power  dependence 


- 61 


Since  the  coefficient  contains  a contribution  from  free 

carriers  which  have  been  calculated  from  a modified  Kane's  model,  we 
might  question  the  accuracy  of  the  carrier  concentration  measurement. 

The  carrier  concentration  was  measured  by  the  Van  de  Pauw  technique, 
which  for  bulk  samples  with  possible  microcracks  and  grain  boundaries  does 
not  have  as  much  accuracy  as  we  would  like.  It  is  not  unreasonable  that  our 
carrier  concentration  measurement  could  be  wrong  by  a factor  of  2.  Since 
c^g  and  c^g  in  Eq.  (4.17)  are  proportional  to  N from  the  form  of  the 
equation  a factor  of  2 error  in  N gives  a factor  of  3 error  in  c^. 

One  car  easily  attribute  part  of  the  error  to  the  experimental  measurement 
process  since  our  crystal  must  have  three  angles  9 , 0 and  p aligned 
correctly.  Moreover,  even  though  we  tried  to  perform  this  measurement  with 
the  laser  unmodelocked  (and  with  no  mode  beating),  fluctuations  in  the  laser 
appear  at  least  three  times  as  large  in  the  third  harmonic  fluctuations. 

The  d^  coefficient  of  CdGeAs0,  as  discussed  in  Appendix  II,  is  accurate 
to  only  ± 40$  . The  ratio  c /d^  is  accurate  to  ± 25$  . 

Since  the  Type  II  THG  tuning  curve  shown  in  Fig.  21  compares  well  with 
theoretical  predictions  for  the  halfwidth,  it  appears  that  we  were  not  thermally 
heating  the  sample.  The  angular  difference  of  12°  (internal  angle)  between 
Type  II  SHG  and  Type  II  THG,  measured  in  the  same  material,  is  in  good 
agreement  with  the  predicted  angle  difference  in  Figs.  3*6  & 3-7  • Consequently 
it  does  not  appear  that  we  were  producing  additonal  free  carriers  to  break 
phasematching.  To  change  the  birefringence  by  0.0024  we  only  need  to  inject 
an  additional  N = 5 x 10  carriers/cnr',  which  would  not  affect  the  value 
of  c as  much  as  indicated  in  our  measurements. 


►V-  — 


m 


By  far  the  most  important  point  to  consider  is  that  our  modified 

Kane  model  is  too  simplified  a picture.  Several  workers  already  observed 

that  Kane | s model  leads  to  predicted  values  of  effective  masses  and  third 

order  susceptibilities  2 or  3 times  smaller  than  actually  observed.  As 

ll8 

noted  in  Appendix  II,  Cardona  in  reflectivity  and  Faraday  rotation 

119 

measurements,  and  Wynne  in  third  order  susceptibility  mixing  experiments, 
both  observed  such  deviations  of  Kane's  model's  predictions  from  experiments. 
In  addition,  Aspnes  has  shown  that  contributions  to  the  second  order 
susceptibility  from  several  points  in  the  Briliouin  zone  must  be  considered. 
He  discussed  the  enhancement  of  the  second  order  susceptibility  from  small 
gap  regions  throughout  the  Briliouin  zone,  rather  than  only  at  the  direct 
gap.  Much  more  work  on  the  band  structure  of  CdGeAs2  needs  to  be  done  before 
these  points  can  be  clarified. 

^ • Efficient  Harmonic  Generation  in  CdGeAs„ 

At  present  all  measurements  in  CdGeAs2  have  been  made  with  200  nsec 
pulses  from  TEA-C02  lasers.  Our  1 nsec  system  described  in  Chapter  II  was 
not  operating  well  enough  to  use  for  measurements.  From  Table  3 we  can 
derive  expressions  for  the  conversion  efficiency  into  the  nth  harmonic  by 
taking  the  ratio  I (noo)/l (cjo) . 

*<2“)  = )(l6^)("~“)2(ti  V2  e*Pt-  + a^ji] 

«?(3cu)  = (-)  (— ) (l6rr  ^ 1 l2(^)exp[-  ■|(a3  + 30^)1] 


- (2)  {')  (lojt  ^ (t£  t^)  I^(cL')exp[-  |(o ^ + ta1)i] 


(^.6) 


Here  tF(ftx))  is  the  power  conversion  efficiency  into  the  nth  harmonic 
and  i(rnj)  is  the  intensity  of  the  nth  harmonic.  In  our  experiments, 

for  example,  we  have  measured  conversion  efficiencies  of  2^  for  SHG  and 

-4 

10  $ for  THG  in  a 2 mm  long  non-AR  coated,  p-type  crystal  with  input 

2 5 

intensities  of  20  MW/cm  . Kildal^  has  measured  a 7$  SHG  conversion  efficiency 

2 

for  a similar  9 long  p-type  crystal  and  no  MW/cm  input  intensity. 

Research  on  CdGeAs_  crystal  growth  at  the  Center  for  Materials  Research 

d 

at  Stanford  indicates  that  this  material  can  be  made  to  grow  n-type  by 
doping  with  indium  or  vanadium.  Moreover,  crystal  growers  at  both  Stanford 
and  Lincoln  Laboratories  have  produced  1 cm  long  pieces  of  CdGeAsg.  Seeded 
vertical  Bridgeman  growth  using  [111]  oriented  seeds  looks  promising  as  a 
technique  for  growing  large  single  crystals.  Also,  by  using  1 nsec  pulses 

r* 

and  focusing  to  intensities  of  100  MW/cm^,  it  is  possible  to  reach  60£  SHG 

and  30=£  THG  conversion  efficiencies  in  1 cm  long  crystals.  As  mentioned  in 

Chapter  I,  burn  densities  in  p-type  Ge  for  a 1.7  nsec  pulse  have  been 

o 

measured  at  1 GW/cni".  From  the  conversion  efficiencies  quoted  above,  if 

2 

the  burn  density  of  CdGeAs^  is  close  to  1 GW/cm  one  would  expect  50  - 80$ 

conversion  to  second  and  third  harmonic  in  a 1 cm  long  crystal  with  tight 

5 

focusing.  However,  Kildal  has  shown  that  there  is  saturation  of  the  SHG 
signal  at  high  power  densities  with  200  nsec  pulses.  The  origin  of 
saturation  appears  to  depend  on  some  nonlinear  absorption  mechanism  involving 
the  SHG  signal.  No  free  carriers  produced  by  multiphoton  absorption  were 
detected;  this  is  the  primary  saturation  mechanism  in  Te,  which  has 

a smaller  bandgap  than  CdGeAs0.  Further  work  on  the  nature  of  the  saturation 
mechanism  in  CdGeAs^  is  essential,  particularly  with  1 nsec  pulses,  in  order 
to  fully  evaluate  its  usefulness  for  high  power  infrared  optics.  We  have 

- 6k  - 


■.■-..•a-i.i.u f 


examined  phasematched  THG  in  liquids  and  gases  as  an  alternative  to  THG 

, 123 

in  CdGeASg,  but  no  experiments  have  been  performed  to  date. 


C.  FOURTH  HARMONIC  THEORETICAL  PREDICTIONS 

In  Section  II  of  Appendix  IV  the  basic  approach  of  the  Bond  Charge 

102 

Model  is  outlined.  Starting  from  Penn's  model  for  the  dielectric  constant 

10?  104 

of  a semiconductor  and  Phillip ' s J and  Van  Vechten ' s theory  for  the 

dielectric  properties  of  semiconductors,  Levine"^"  developed  the  Bond  Charge 

Model  to  calculate  second  order  nonlinear  susceptibilities  of  binary  semi- 

105 

conductors.  Chemla  later  extended  this  model  to  ternary  compounds,  and 
the  paper  given  in  Appendix  II  is  a generalization  to  third  order  suscept- 
ibilities. 


In  this  section  we  apply  the  Bond  Charge  Model  to  calculate  the 
fourth  harmonic  tensor  coefficients  in  CdGeAs^.  We  begin  by  writing  the 

polarizability  of  a bond  in  the  language  of  Phillip's  mean  energy  gap  E 


g 


2 

o D E 

a . - (a(|  + 20^)  . (2ac)  — § 


(4.7) 


2 2 —9  4 2 

where  aQ  = h /me  = 5*3  X 10  cm  , Eq  = me  /2  h = I3.6  ev.  D is  a 

correction  factor  near  unity  to  account  for  effects  of  core  electrons  and 


E is  given  by 
8 


E2  = H2  + C2 
g 


(4.8) 


The  honopolar  and  heteropolar  energy  gaps  H and  C are  given  in  Fnillip's 


theory  by 


2 

CAB  = 1'5e  «P 

rA  rB 


; I VrA  + rB>] 


X “ (rA  - rCA,?  + (rB  ' rCB)5 
AB 


(4.9) 


where  zft  , rA  , r^  are  the  valence  number,  covalent  radius106  and  the 

core  radius  of  the  A atom  in  the  AB  bond,  and  k is  a Thomas-Fermi 

s 

screening  constant.  Higher  order  susceptibilities  arise  from  nonlinear 
variations  of  the  mean  bond  polarizability  a with  respect  to  electric 
field.  We  can  write  the  dipole  of  a bond  in  terms  of  a Taylor  series 
expansion  of  a ; 


Pi  - “il  Ei 


?i  (ai^E  = 0+(bE)lE  = 0Ej  + ””  Ei 


(4. 1C) 


The  derivative  terms  give  the  nonlinear  susceptibility  components  for  each 
bond.  For  the  case  of  FKG  we  have 


Jk'P  3 1 iEj  E = 0 


(4.11) 


As  discussed  in  Appendix  II,  we  assume  that  we  can  talk  about  these  higher 
order  polarizabilities  in  terms  of  a bond  charge101  q , located  at  distances 
rA  anl  rB  ^rom  each  atom  in  the  bond.  We  calculate  the  variations  in  a 
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with  electric  field  by  examining  the  variations  in  rA  and  rg  with 
electric  field.  It  is  shown  in  Appendix  II  that 


Ar  = 

A 


ArA  = 


-SrB  . 

B g 


<££. 


4* 


— 5 ; 4i'b=  ttj 

2rAq  2rBq 


E parallel  to  bond  axis; 


E perpendicular  to  bond  axis. 


(4.12) 


Since  the  variations  in  rA  and  rB  are  second  order  in  E^  , we 

retain  only  even  powers  of  E.  . Consequently  the  odd  order  derivatives  of 

X 

with  respect  to  Ei  vanish: 


2p  + 1 


(4.13) 


All  order  derivatives  in  Ej  | are  retained.  By  using  Eq.  (4.13)  and 
Kleinman  symmetry  one  can  show  that  the  elements  of  the  fourth  order 
susceptibility  of  a bond  are  given  by 


6 = 6 
zzzzz 


1 

b E|  I 


1 6 «L 

= 6,  = ~ — = 0 


xx xxx  I 6 2 

X y E. 


3 

5 6 1 u XZ 

xxzzz  = I II  ~ 6 XE  6e  f 


+ all  permutations 


(4.14) 
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By  combining  expressions  Eq.  (4.7)  through  (4.14)  as  discussed  for 
the  case  of  THG  in  Appendix  II,  we  obtain  an  expansion  for  ^ j (AB)  for 
the  fourth  order  polarizability  of  the  bond  between  atoms  A and  B 
given  by  the  following: 


6, i (AB)  = 


-30=' 


AE 


4(2fh  - Dh^H  + 4(2f.  - l)gjc  + 4(6f.  - l)hlg^H 

+4(6fn  - l)h^g]C  - 8(h1g2  + h^jHC 

-£^2(4^  - l)h1h2  + 2(4f.  - l)glg2  - 2(y  fg  C)] 


(4.15) 


where 
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^2  = ~2~ 
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=rA  + rB,‘ 
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. 2 
-be 
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— r — r 

exD  1"  — ksfr  + r ^ 1 

-103.2 

h 
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3 q3 
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Here  q is  the  bond  charge  given  by 

2 

_ i 

q 


■-3  [4 


+ 0.4(1  - f. ) e 


i>} 


(4.16) 
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and  the  factor 


where  n is  the  index  of  refraction,  f^  the  ionicity^^ 

3 arises  from  the  fact  that  ve  have  a tensor  derivation.  In  CdGeAs^  we 
must  derive  the  value  of  the  polarizability  for  the  CdAs  bond  and  the 
GeAs  bond,  denoted  by  6^  and  6^  respectively. 

Finally  we  must  sum  over  the  contributions  of  all  bonds  in  the  unit 
cell  of  volume  V.From  Fig. 3-1  we  have  two  nonzero  components  of  the  macro- 
scopic FHG  tensor.  In  terms  of  the  parameters  t and  a , which  are  a measure 
of  the  distortion  of  the  unit  cell,  we  can  express  b^^  and  b^,_^  as 
the  geometrical  sum  over  the  bond  polarizabilities 


16 


3K4) 


9\/3V 


hy-f 


\ + /,  , 4t  2a  \ - 

) 6n  - (!  3 + 3 ) 5 n 


16 


*1(5) 


9 J3V 
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(1  . 2t  a}  6+ 

v 3 3 ' 11 


i] 


(4.17) 


_ c kx 

where  t = — a = ~ — 1 , and  a and  c are  the  lattice  constants.  By 

substituting  the  expressions  Eq.  (4.15J  for  and  into  Eq.  (4.17) 

(using  values  of  the  parameters  for  CdAs  and  GeAs  bonds  found  by  interpolation 

IOI.I03-IO6, 


of  values  for  other  bonds," 


')we  derive  the  values 


| bl (4 ) | = 2.80  x 10- 18  esu 
l^l(5)l  = O.72  X 10  esu 


Unit  s cm 


/ 


/ erg3/2 


(4.18) 


This  is  the  central  result  of  this  section.  The  only  experimental  results 

-21 

on  phasematched  FGH  in  any  material  give  a value  of  b « 10  esu  for 

124 

lithium  formate.  The  large  nonlinear  response  of  CdGeAs^  is  apparent 

in  the  three  orders  of  magnitude  difference  between  the  two  susceptibilities. 

Finally  we  can  estimate  the  conversion  efficiencies  to  be  expected  for 
8 2 

FHG  in  CdGeAs^  with  10  W/cm  input  intensity  in  a 1 nsec  pulse.  We  assume  a 
1 cm  long,  AR  coated  crystal.  From  Eq.  (i+.6)  we  calculate  a conversion 

-8  i 

efficiency  of  10  <&  for  a typical  p-type  crystal  (a^  =0.5  cm'1  , 

\ = 7-0  cm'1). 

For  a high  resistivity  sample  (c^  ^ 0.01  cm'1  , « 4.6  cm'1),  and  an 

input  intensity  of  lGW/cm  the  conversion  efficiency  should  be  about  10~^ 
for  FHG.  It  appears  that  FHG  in  CdGeAs^  will  never  be  a practical  source 
of  infrared  radiation,  but  coupled  with  SKG  and  TKG  results,  it  may  provide 
additional  insight  into  the  theoretical  models  for  nonlinear  susceptibilities. 
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APPENDIX  I 


CALCULATION  OF  NONZERO  ELEMENTS  OF  nth  RANK  SUSCEPTIBILITY  TENSORS 

In  this  appendix  we  examine  a very  straightforward  procedure  for 
determining;  the  nonzero  elements  of  a nonlinear  optics  susceptibility 
tensor  to  any  order  from  one  of  lower  rank.  Since  most  work  in  the  past 
has  required  consideration  of  only  , x/2^  and  sometimes  x^)  tensors, 

the  direct  inspection  method  ■ has  been  most  often  used  to  simplify  the 
forms  of  tensors.  However,  one  can  easily  see  that  this  method  can  become 
quite  cumbersome  for  higher  rank  tensors.  The  total  number  of  elements  in 
a tensor  rank  n is  3°  , corresponding  to  a susceptibility  of  order 
(n-l),  namely  . For  fourth  harmonic  generation  for  example,  n = 5 

and  3 = 2^3  elements.  For  a possible  fifth  order  process  n = 6 and 

3^  = 792  elements.  Nevertheless,  the  direct  inspection  method  has  been 

3 

applied  to  fifth  and  sixth  order  tensors. 

From  another  point  of  view,  one  can  apply  the  methods  of  group  theory 
to  the  task  of  expressing  a tensor  of  rank  n in  terms  of  its  irreducible 
components  of  lower  ranks,  1 This  technique  has  been  used  successfully 

in  nonlinear  optics  by  Jerphag.ncn.  The  use  of  such  techniques,  particularly 
the  realization  that  one  can  very  closely  associate  the  properties  of  higher 
rank  cartesian  tensors  with  higher  order  spherical  harmonics,  leads  to  a 
very  elegant,  but  also  very  practical  formulation  of  higher  order  tensors 
for  many  physical  problems.  - This  method  moreover,  provides  the  only  simDle 
way  for  finding  the  exact  relationships  between  the  nonzero  elements  of  a 
tensor  of  rank  n . 
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In  this  appendix  we  simply  borrow  some  points  from  these  more 
general  treatments  and  consider  the  case  of  nonlinear  optical  susceptibility 
tensors.  There  are  three  classes  of  symmetry  conditions  we  must  consider: 
the  symmetry  of  the  particular  crystal  class  in  question;  overall  permutation 
symmetry;1^  and  Kleinmen  symmetry.15^  For  a tensor  of  rank  n and  a 
harmonic  generation  process,  n - 1 of  the  frequencies  involved  are  the 
same.  Overall  permutation  symmetry  allows  one  to  permute  n - 1 of  the 
tensor  element  subscripts  without  changing  the  value  of  the  coefficient. 


For  SHG  for  example 

dijk  " dikj 


(1.1) 


If  all  of  the  frequencies  in  question  are  within  the  transparency  range 
of  the  material,  one  often  can  neglect  dispersion  in  the  nonlinear 
susceptibility.  In  that  case  Kleinman  symmetry  holds,  and  one  can  permute 
all  n subscripts.  Again  for  SHG  we  have 


For  a tensor  of  rank  n , the  application  of  Kleinman  symmetry  makes  the 
tensor  fully  symmetric.  The  total  number  of  independent  elements  is  given 


For  overall  permutation  symmetry  the  tensor  is  symmetric  in  only  n - 1 


1 

•1 

: 

f 


indices.  The  total  number  of  independent  elements  is  given  by 


127 


3 (VA)  - 

Finally,  from  above,  for  the  case  of  no  special  symmetries  there  are 
independent  elements.  The  factor  3 , of  course,  comes  from  dealing  with 
a space  with  three  dimensions  x,  y,  z.  Table  1.1  summarizes  the  number 
of  independent  elements  for  different  symmetry  conditions,  before  the 
application  of  the  properties  of  the  crystal  class.  Kleinman  symmetry  is 
a very  powerful  condition  for  reducing  the  number  of  independent  tensor 
elements  and  it  should  always  be  applied  first  before  the  symmetry  properties 
of  the  crystal  are  considered. 

Next,  we  can  write  a tensor  of  rank  n in  two  ways:1^  first  it  is 
the  tensor  product  of  n vectors 


(n 

ijkl P 


© d1'® 


(I.?) 


Secondly,  it  can  be  written  as  the  tensor  product  of  a vector  with  a 
tensor  of  rank  (n  - 1) 


(n 

ijkl P 


c(n_-0 

jfcl p 


(1-4) 


In  particular,  one  can  write  each  higher  order  nonlinear  optics  tensor  as 
the  tensor  product  of  a vector  and  the  previous  lower  order  tensor. 


X. 


(?) 

ijk 


h ® x3k’ 


X,. 


(3) 

ijkl 


h ® xju 


(I.! 
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The  prescription  of  ‘■his  0pp**ndix,  therefore,  is  to  first  apply 
Eq.  (1.4)  to  generate  the  appropriate  nonzero  elements  of  the  n rank 
tensor,  but  using  the  Kleinm^n  symmetry  reduced  form  of  the  (n-1)  rank 
tensor.  Then,  the  properties  of  the  particular  crystal  symmetry  can  he 
applied  to  the  already  much  reduced  tensor.  As  a practical  example,  we 
can  generate  the  independent  elements  of  y[^  and  yj3^  . We  first 
note  that  the  linear  susceptibility  has  six  independent  elements 


(11  22  33  3?  13  12) 


,(1) 


(1.6) 


where  1,  2,  3 ore  short  notati'  for  x,  y,  z.  To  generate  y(?'  we 

ijk 

write 


ijk 


v.(S>  - t'1)  ® 


(123)  X (11  22  33  32  13  12) 


(1-7) 

If  we  write  down  all  the  possible  combinations  in  the  product  and  keep 
only  the  independent  ones  we  have 


(111  112  II3  122  123  I33  222  223  233  333) 


:(2) 

ijk 


(1.8) 


Note,  there  are  ten  independent  coefficients  as  predicted  in  Table  1. 1 
for  Kleinman  symmetry.  Using  Eq.  (1.8)  we  can  determine  the  coefficient 


of  *(3)  . 

ijki 


x(3) 

ijki 


* t(l)  © XC2) 

i w jk£ 


(1-9) 
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TABLE  1.1 


INDEPENDENT  ELEMENTS  OF  A TENSOR  OF  RANK  n 


Tensor 
Rank  (n) 


Susceptibility 
Order  (n-1) 

No  Symmetry 

(3°) 

1 (to) 

9 

2 (2co) 

27 

3 (3o)) 

81 

u (M 

2U3 

5 (5oi) 

729 

Overall  Permutation 
Symmetry  ■3n_(n'*'1 ) 


Kleinma.n 
Symmetry 
/ (n+l)  'n+°)\ 

\ 7 ) 


2 

3 

k 

5 

6 


18 

30 

U5 

63 


6 

10 

15 

21 

23 


Whiting  out  the  product  and  keeping  the  independent  elements  gives 


I 


Here  there  are  fifteen  independent  coefficients,  again  in  agreement  i-ith 
Table  1.1.  At  this  point  we  could  apply  the  crystal  class  symmetry 
properties,  i.e.  h2  m symmetry  fo~  CdGeAs2  and  obtain  the  final  expressions 
for  vijk  and  vijki  " We  would  obtain  the  tensors  shewn  in  Fig.  10. 

The  procedure  outlined  here  wee  used  to  calculate  the  FHG  tensor  for  CdGeAs2 
also  shown  in  Fig.  10.  The  alternative  approach  of  starting  with  the  21*3 
element  general  tensor  and  using  the  direct  inspection  method  firsv,  oroved 
too  awesome  to  this  author.  Finally,  this  technique  was  used  to  calculate  the 

nonzero  independent  elements  of  the  fifth  harmonic  tensor,  although  it  is  rot 
shown  in  this  manuscript. 


A3TCNDIX  II 


EXIERIMEUTAL  AND  TIEORETICAL  STUDIES  OF  THIRD  HAH'CNIC  GENERATION  HI 

THE  CHALCOPfRITE  CdGeAs2 

The  extraordinarily  large  third  order  susceptibility  measured  for 
CdGeASg  in  THG  experiments  provoked  a detailed  investigation  of  the  origins 
of  higher  order  nonlinear  susceptibilities  in  semiconductor  crystals.  The 
paper  provided  in  this  appendix  summarizes  the  essential  features,  both 
experimental  and  theoretical,  of  the  large  third  order  response  in  CdGeAs^* 
The  three  principal  results  of  this  paper  are  thp  following: 

1.  All  nonlinear  optical  susceptibilities  are  due  to  enharmonic 
response  of  electrons.  However,  by  symmetry  arguments  one  can  conclude 
that  even  order  susceptibilities  should  reflect  properties  dep»  nding  on 
crystal  symmetry  alone.  Consequently  they  should  only  involve  electrons 
localized  in  the  crystal  bonds.  Odd  order  susceptibilities,  on  the  other 
hand,  are  still  present  In  isotropic  media.  Consequently,  the  isotropic 
free  carrier  cloud  should  exhibit  odd  order  response,  whereas  its  contrib- 
ution to  even  order  susceptibilities  should  be  zero.  This  suggests  the 
interesti..<j  possibility  of  providing  direct  control  over  the  magnitude  of 
higher  order  susceptibilities  by  controlling  carrier  concentration.  For 
even  order  susceptibilities  we  wish  to  hove  a high  resistivity  material 
to  eliminate  as  many  absorption  mechanisms  as  possible.  For  odd  order 
susceptibilities  the  situation  is  not  so  obvious.  Larger  free  carrier 
concentrations  lead  not  only  to  larger  nonlinear  coefficients,  but  usually 
to  larger  ebsorption  losses  as  well.  At  seme  concentration  there  is  an 
optimum  nonlinear  response. 
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2.  CdCteAs  h*s  ep.ouch  birefringence  to  allow  both  Type  I and  Type  II 
THG.  A careful  examination  of  nonzero  elements  of  third  harmonic  tensers 
for  various  symmetry  materials  reveals,  discussed  in  Section  IV,  that 
one  can  separate  bond  electron  fro r.  free  carrier  effects  by  phascmstchinc 
experiments  alone. 

3-  The  measurement  of  several  higher  order  nonlinear  optical 
susceptibilities  in  the  s^re  material  provides  a new  technique  for  studying 
some  properties  of  electron  charge  clouds  and  bend  structure  in  semi- 
conductors. All  order  susceptibilities  have  contributions  from  bond 
electrons,  which  are  proportional  to  higher  order  multipole  moments. 

Such  an  investigation  could  provide  further  insight  to  the  models  currently 
used  to  describe  nonlinear  optical  response,  such  ps  the  Bond  Charge  Mod'_i 
or  the  Bond  Orbital  Model.  Free  carriers,  on  the  other  hand,  are  responsive 
to  nonparebolic  energy  bands  and  provide  a means  cf  testing  Kane's  Model 
for  energy  band  shapes,  for  example.  These  points  have  been  discussed 
further  in  Section  IV  and  are  alluded  to  in  the  following  paper. 
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Experimental  and  Theoretical  Studies  of 
Third-Harmonic  Generation  in 
the  Chalcopyrite  CdGeA^ 

DANIEL  S.  CHEMLA.  R.  F.  BEGLEY,  and  ROBERT  L.  BYER 


Abitmcl-  Experiiaewtit  tud  theoretical  studies  of  rhM-Wmuriic 
(me ration  (THC)  in  the  cbalcopvrile  semiconductor  CdGeAs,  arc 
presented.  The  ptasc-mitchinf  confip.itions  for  THC  arc  anal  vied  from 
the  irreducible  components  point  of  srw.  A thenrv  of  the  bound  rlfctroo  and 
frro-carricr  contribution  lo  the  third-*  <ter  ui-.-rotsility  it  presented.  Th» 
experimental  results  arc  given.  The  effective  nonlinear  -oelheient  for  tvpe-ll 
THC  It  mainly  due  'j  the  frrr -carrier  contribution,  and  for  a hole  roncen- 
tratlon  of  S x tOu  cm'  it  is  measured  to  be  H3  - 6)10  ■ LSI . The  prac- 
tical applications  of  THC  in  CdGeAs,  are  divussed. 

I.  Introduction 

SINCE  the  first  den onstration  of  phase-matched 
second-harmonic  gereration  (SHG)  in  AgGaS,  [I] 
the  second-order  optical  effects  of  the  ternary  compounds 
with  chalcopyrite  structure  have  been  widely  studied 
(2)-(6J  Their  large  nonlinear  susceptibilities  together  with 
their  natural  birefringence  make  them  very  attractive  for 
nonlinear  optical  devices.  Nonlinear  optical  mixing  has 
been  recently  achieved  in  ZnGeP,  (7).  AgGaS,  [8],  [9],  and 
AgGaSe,  [10]  High  efficiency  second-  ?nd  third-harmonic 
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R 1 Beglev  and  R.  L Bvcr  are  with  the  Microwave  Laboratory.  W. 
Hansen  Laboratoties.  Stanford  I'mversny.  Stanford.  Calif  94)0? 


generation  (THG)  has  been  demonstrated  in  CdGeAs, 
[II],  [12J. 

In  addition  to  its  large  second-order  nonlinearity 
CdGeAs,  has  a correspondingly  large  third-order  suscep- 
tibility. In  this  p?pcr  we  report  both  theoretical  and  ex- 
perimental studies  of  the  third-order  susceptibility  of  the 
II-IV-V,  chalcopyrite  compound  CdGeAs,. 

The  third-order  susceptibility  is  a tensor  of  rank  4 
(1 3]— [1 5J  In  the  case  ofTHG  this  tensor  <•.,»,  is  symmetric 
in  the  three  last  indices  (16)  and  a contracted  notation  can 
be  used.  Thus  c,,„  — c,m  wherem  runs  from  0to9  |I6]  and 
the  nonlinear  polarization  can  be  written  as 

fi  * c.tii(3u),  w,  u,  u)E,'’ Et" Ei“  ■ (I) 

The  Einstein  summation  convention  is  assumed  and  8m  is  a 
ten-element  row  matrix  defined  in  Appendix  I 

In  a crystal,  far  from  the  electronic  and  vibrational  ab- 
sorption regions,  the  susceptibilities  exhibit  a very  small 
frequency  dependence.  Therefore,  the  Kleinman  [17]  con- 
jecture is  valid  and  c„»,  is  independent  of  the  index  per- 
mutation so  that  only  15  of  the  81  components  are  in- 
dependent. 

The  intrinsic  svmmetrv  of  a tensor  is  better  understood 
when  use  is  made  of  its  decomposition  into  a sum  of 
irreducible  components  (IS).  The  decomposition  of  the 
symmetric  fourth-rank  tensor  c„„  involves  [19]  a scalar 
(weight  0.  one  component)  a deviator  (weight  2.  five  com- 
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ponents),  and  a nonor  (weight  4,  nine  components)  The 
scalar  ?*  - £„  <•„„  is  nonzero  for  all  media  and  describes 
the  geometrical  average  of  the  Cartesian  components  The 
deviator  c • £,  c(/„  - is  nonzero  for  the 

uniaxial  and  biaxial  crystals  and  describes  a difference  of 
behavior  between  the  crystal  axes,  analogous  to  the 
birefringence  The  nonor  - r„„  - 1/7  I.XI, 
c„*<)  + I '3V4.4«(  ♦ + b„l, »)(£,.,  <■"»•)  describes  the 

departure  from  isotropy  In  the  expression  of  c„t,* 
denotes  the  summation  ever  the  distinct  permutations  of 
the  four  indices  The  three  irreducible  components  obey 
different  transformation  laws  and  therefore  have  diffc'ent 
contributions  to  the  effective  third-order  susceptibility  un- 
der phase-matching  conditions 
Because  the  point  group  of  CdGeAs,  is  42m.  the  third- 
order  susceptibility  tensor  can  be  written  in  the  contracted 
notation  [20] 


Cll C|fC„  • •"* 

C|l  Cli  • • • • C,»- 

.*  ■ C||-C|#-C|A  •••  . 


(2) 


According  to  the  Midwinter  and  Warner  classification  [20] 
there  can  be  three  types  of  phase-matched  third-harmonic 
generation  (PM-THG)  in  a crystal  belonging  to  the  class 
42m.  The  corresponding  effective  nonlinear  coefficients  for 
a positive  birefringence  crystal  are  [20] 


fi  " "Ifru  3c,«)  cos  ( sin  4$,  e -f-  t -f-  e — > o 

(3a) 

Cii  “ Il(Cii  — 3c, ,)  cos’  #sin’  2# 

+ (c,»  — c„)  cos’  6 + c,tJ  sin1  a cos  a. 

e + e -f  o — * o (3b) 


Cm  " [}(c„  — 3c„)  cos  8 sin  4$]  sin  a cos’  a, 

c -f-  o + o — ♦ o (3;) 


where  the  angles  6 and  d define  the  direction  of  propaga- 
tion for  phase  matching  in  the  ci  vstalographic  axes  X YZ. 
and  a is  the  angle  between  the  fundamental  field  and  the 
XY  plane. 

The  combination  (c„  - 3c„)  is  a component  of  the  no- 
nor. It  characterizes  the  departure  from  isotropy  in  the  X Y 
plane  The  combination  (r„  - <*,«>  describes  the  difference 
between  the  direction  Y and  Z and  appears  in  the  deviator. 

It  is  well  known  from  the  theoretical  [13],  [2I]-[26]  and 
experimental  [27]-[30]  studies  of  third-order  nonlinear 
processes  in  elementary  and  binary  compounds  that  in 
semiconductors  the  third-order  susceptibility  arises  from 
two  distinct  contributions 


Ci. 


+ o-* 


(4) 


where  f,*rc  is  the  contribution  from  the  free  carriers  and 
c,m*f  is  the  contribution  from  ihe  bound  electrons 

In  the  crystal  the  bound  electrons  hjve  the  symmetry 
42m  and  contribute  to  all  the  effective  nonlinear 
coefficients  However,  the  free-carrier  gas  has  a higher 
symmetry.  For  example,  the  free-carrier  contribution  to 
the  second-order  susceptibility  vanishes  in  the  dipole  ap- 
proximation even  in  a noncentric  semiconductor  In  a 
crystal  w ith  point  group  42/m  the  tree  carriers  are  expected 
toheisotropicintheA  ) planeandthereforer,,'*  • V,/r  In 
type  I and  type  III  PM-THG  only  the  bound  electronscon- 
tribute  to  the  effective  nonlinear  coefficients,  whereasin  type 
II  PM-THG  both  the  bound  electrons  and  the  free  carriers 
contribute  It  snould  be  noted  that  this  result  isquitegeneral 
because  in  the  maioritv  of  crystal  point  groups  the  free- 
carrier  gas  is  isotropic  or  nearly  isotropic  in  the  X Y plane 
The  requirement  of  invariance  under  axial  symmetry  on  a 
tensor  of  rank  n isequivalem  to  the  invariance  under  a rota- 
tion 2ir/m  w ith  mi  >n  In  thecaseofTHG  the  relevant  tensor 
is  of  rank  4.  The  possible  phase-matching  configurations  for 
the  fret -carrier  gas  are  the  same  as  for  say.  an  hexagonal 
crystal.  From  [20,Tablc  I ] it  can  b.  seen  that  in  thiscasetvpe 
I and  type  III  PM-THG  are  impossible.  In  the  majority  of 
semiconductorsthe  free-carrier  contributions  PM-THG  is 
of  type  II. 

In  all  the  previously  reported  determination  of  the 
third-order  susceptibilities  of  semiconductors  onlv  mixing 
experiments  (w,  + u,  - w,  ■ w,  with  w>  » w,|  had  a sub- 
stantial coherence  length  Thus  scleciive  phase  matching 
could  not  be  used  to  separate  the  two  contributions  to  cim. 
The  two  effects  were  separated  by  varying  the  carrier  con- 
centration. In  Cri' je.As,  both  type  I (e  e + t — o)  and 
type  II  (e  -»■ t + o — o)  PM-THG  are  abowed  so  that  free 
carrier  and  bound  electron  contribution  to  c,m  can  be 
separated  experimentally  for  the  first  time  by  choice  of 
crystal  orientation 

Another  aspect  of  this  s'mmetrv  analysis  is  that  it 
provides  a useful  guide  to  evaluate  the  relative  magnitude 
of  some  irreducible  components  of  unknow , high-rank 
tensors  like  c,m.  The  i^edueible  components  1 ->wer 
weight  of  the  tensor  cim  are  a scalar  and  a deviatu  The 
bound  election  contribution  is  related  to  the  po'arizaoilitv 
[22J  of  the  bonds  and  the  tree-carrier  contribution  is 
related  to  the  derivatives  of  the  energy  with  respect  tr  the 
wavevector  [13].  From  the  linear  optical  and  transpa. ' 
properties  of  the  crystal  it  is  eass  to  form  a scalar  and  a 
deviator,  respectively,  related  to  the  bond  polarizability 
and  the  derivative  of  the  energy  Therefore,  an  evaluation 
of  the  ratio  (r„  - c, ,)'<•, t can  be  given  for  the  two  systems 
of  charges  within  the  crystal  For  the  bound  electrons  v,c 
expect  the  ratio  to  be  of  the  order  of  magnitude  of  (*  - 
Xi)  X«*  where  x and  are  the  linear  susceptibility  com- 
ponents. For  the  free  carriers  this  ratio  should  be  of  the 
same  order  of  magnitude  I's  (m  — /n.Vni,  wherem  and 
nii  are  the  effective  masses  of  the  carrie-s  For  a p type 
CdGeAs,  the  ratio  is  of  the  order  of  unity  for  an  n-ivpe 
material  it  is  of  the  order  of  0.5  according  to  the  elfecttve 
masses  given  in  [31] 
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The  contribution  of  the  bound  electrons  is  examined  in 
Section  II  In  that  section  we  generalize  the  'bond  charge 
model”  (32J  to  a tensor  form  which  enables  us  to  calculate 
•I!  the  components  of  the  third-order  polarizability  of  a 
tetrahedral  bond  with  axial  symmetry 
In  Section  III  we  calculate  the  contribution  of  the  free 
carriers  This  is  performed  by  using  Kildal’s  (31 1 
generalization  of  Kane's  model  (33]  to  crystals  with 
chalcopynte  structure 

In  Section  IV  the  experimental  work  is  presented.  The 
effective  nonlinear  coefficients  corresponding  to  type-1  and 
type-ll  PH-THG  have  been  compared,  cn  has  been 
measured  relative  to  </u*  by  comparing  PH-THG  to  the 
mixing  of  the  laser  light  with  tvpe-ll  PM-SHG  in  the  same 
crystal  but  at  different  orientation  i#,,**0  - tf„™°  - II* 
± I*)  and  to  d, 4 by  comparing  PM-THG  to  PM-SHG 
Finally,  the  comparison  of  experimental  and  theoretical 
results  are  discussed  The  major  practical  consequence  of 
our  work  is  to  point  out  the  possibilitv  of  high  conversion- 
efficiency  type-ll  PM-THG  in  heavily  doped  chalcopynte 
semiconductors. 

II  Theory  of  the  Third-Order  ScscErnaiLmr  of 
the  Bound  Electrons  tN  a Tetrahedral  Crystal 

The  most  successful  recent  theory  of  the  dielectric 
properties  of  elementary  and  binary  tetrahedral  solids  is 
due  to  Phillips  (34]  and  Van  Vechten  (35],  Their  theory 
has  successfully  been  extended  to  multibond  tetrahedral 
solids  (36]  and  is  the  basis  of  several  models  describing  the 
second-order  nonlinear  susceptibilities  of  various  crystals 
(32].  (37]— (40].  We  present  in  this  section  a generalization 
of  the  bond-charge  model  (32]  to  calculate  all  the  com- 
ponents of  the  nonlinear  optical  susceptibilities  of  difTc’ent 
orders,  and  we  consider  in  more  detail  the  case  of  the  third- 
order  susceptibility. 

To  give  a good  description  of  all  the  components  of  the 
susceptibilities  we  have  to  keep  track  of  the  tensor  form  of 
the  different  quantities  involved  in  the  theory  The  bond  in- 
terpretation of  the  Phillips  and  Van  Vechten  theory  (36) 
shows  that  a mean  bond  polarizability  can  oc 
macroscopically  defined  in  a way  • hich  includes  the  local 
field  effects  by 


r * 

or  » ](<*,  + 2ox)  w (2flo)  D — y (3) 

In  (5)  a and  c»i  are  the  parallel  and  transverse  com- 
ponents of  the  axially  symmetric  bond  polarizability  ten- 
sor, r,  i$  the  Phillips  mean  energy  gap  (34],  (33] 

£,'  - £.*  + C*  (6) 

D is  a correction  factor  close  to  unitv  which  describes  the 
core  electron  contributions  [V3|.  and  a * ~ me ’ is  the 

Bohr  radius  and  t„  - me*  ?Ji:  is  the  hvdrorcn  ionization 
energv  In  (."■)  the  Penn  correction  factor  has  been 
neglected 
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Fig.  t.  V aha  t ton  of  (He  bomJ-chirft  position  under  an  appticd  field 


Following  Levine  (32]  we  assume  that  the  nonlinear 
susceptibilities  arise  from  the  variations  of  the  linear 
susceptibility  due  to  the  applied  field  t The  microscopic 
dipole  of  the  bond  p,  * o,i(*)e,  can  be  written  by  expand- 
ing the  linear  polarizability  (41) 


This  equation  gives  the  relation  between  the  different  com- 
ponents of  the  bond  higher  order  polarizabilities  and  the 
derivatives  of  the  linear  polarizability  components  with 
respect  to  the  applied  field 


It  is  important  to  note  that  the  derivatives  are  evaluated  at 
vanishing  field  so  that  the  model  is  dispersionless. 
Therefore.  Klcinman  symmetry  (17]  holds.  i,e.. 
y,  _ • • arcinvariantunderpermutationsoftheirindices. 

It  is  now  ncressarv  to  calculate  the  field  dependence  of 
the  a,,.  In  the  tond  charge  scheme  |32]  the  motion  of  the 
weakly  bound  aid  well-localized  bond  charge  (42],  (43]  q 
is  responsible  for  the  bond  polarizabilities.  The 
homopolar  part  Ek  and  the  heteropolar  part  C of  the  mean 
energy  gap  Et  of  the  biatomic  bond  AB  are  functions  of 
the  covalent  radii  [44]  of  the  atoms  rA,  />. 

If  we  interpret  rA  and  r B as  the  distances  between  the 
bond  charges  and  the  two  atoms,  then  the  variations  of  a 
can  be  calculated  through  «he  variations  of  rA  and  Two 
configurations  have  to  be  considered,  as  shown  in  Fig.  I, 
depending  on  the  direction  of  the  field  with  respect  to  the 
bond  axis.  First  the  applied  field  is  parallel  to  the  bond 
axis,  then  A/*  = - Az>  = (o  < ,/q)  are  first  order  in  the 
field.  In  the  second  case  the  field  is  perpendicular  to  the 
bond  axis  and  Ar*  = (ai't'/lrrf)  and  Arfl  ■ 
(a,**.*  arc  second  order  in  the  field  Therefore,  in 

the  expansion  of  the  linear  polarizability  in  powers  of  the 
applied  licid  only  the  even  powers  of  can  appear 
Whatever  is  the  dependence  of  the  au  in  the  field,  their 
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odd-order  derivatives  m «,  vanish 


whereas  the  odd-  and  even-order  derivatives  in  < can  be 
nonzero 

This  remark  together  with  the  Klemman  svmmetry  [17] 
conditions  allows  the  calmlations  of  all  the  components  of 
the  bond  hyperpolanzabilmes  from  the  knowledge  of  the 
mean  linear  polarizability  alone 

Ixt  us  consider  in  detail  the  case  of  the  second-  am' 
third-order  bond  polarizabilities  The  nonzero  corr.- 
oonents  of  the  polar  zabilmes  of  the  axially  symmetric 
bond  are 


itii  km.  »v«t  or  uuantim  iuciiomcs.  mhmi  1974 
TABLF  I 

ExTRHIIOS  Of  THI  Ql  ASflflH  DiHMVi  llll  DirIVxIIIM  Of  TMf 

Howotolaa  asd  HiTmorrt  a«  Pa»to»  tim  Mias  Mono  Emxgt 
Gap  ion  a Hiatomic  Bomj  AB 


i parallel  lo  the  bond  am 


i perpendicular  io  the  bond  am 


fit  m I*...  — fi,„ 

7i  7«*i  7i  “ 7#m  * yy....  7u  " 7mi  • 

Using  (8)  and  the  Kleinman  symmetry  to  rearrange  the  in- 
dices 


(9b) 


From  the  Phillips  and  V*n  Vechten  theory  we  know  a ■ 
(aM  + 2ari)/3  and  therefore  using  (5).  (8).  ind  (9) 


fit  m 3 


• -• 


7. 


(10) 

The  ab  initio-calculations  of  Flytzanis  and  Ducuing  [45] 
have  shown  that  3i/0  « I.  In  the  bond-charge  model 
this  quantity  goes  to  zero.  Ft  should  be  noted  that  in  the 
previous  derivations  of  the  bond  second-order 
polarizabilities  [32]  the  factor  3 arising  from  the  tensor 
form  of  a„  and  0,y»  has  been  omitted.  The  fact  that  the 
linear  polarizability  is  nearly  .spherically  symmetric  is  mis- 
leading: even  if  a =•  a a < aL  their  derivatives  are  com- 
pletely different.  In  the  following  we  will  use  the  relation  a 
* «n  * a*,  but  carefully  evaluate  their  derivatives  ac- 
cording to  (9)  and  (10). 

Let  us  write  the  expansion  of  the  homopolar  and  the 
heteropolar  part  of  the  mean  energy  gap  in  the  following 

way. 


£»(«)  *■  E*  + (a«)A,  + (a<)’A>  + • • • (lla) 


C(0  - C + (or*)*,  + (a*)’*,  + • • • . (Mb) 
The  quantities  h,  and  g,  are  given  in  Table  I for  tne  two 


configurations  considered  Then  expanding  the  mean 
linear  polarizability  by  using  (5)  and  (II)  gives 


+ i Kf.  - I)*,’  + t4 /.  - I)z,’ 

+ 8 (/./.)  '’*,*,  - 2 lh,E„  + t,C))  X a*.’  + • • (12) 

where/,  is  the  ioniciiy  and/,  - I -/.  In  deriving  (12)  we 
have  used  Levine’s  [32J  definition  off,.  and  a is  evaluated 
according  to  (15). 

The  relevant  components  of  the  bonr  polarizabi  x 
are  easily  caiculaleu  from  ( 12),  (10).  and  Table  I.  Equa'  .*n 
(12)  has  been  written  in  a way  which  makes  directlv  ap- 
parent the  ionic,  the  covalent,  and  the  mixed  contributions 
to  the  hyperpolanzabilmes. 

Since  the  microscopic  analogs  of  the  Miller  s delta  [46] 
(i.e.,  0'a'  ano  */a *)  are  very  simple  functions  of 
macroxcopically  defined  quantities  like /,  and  rA  and  rt  the 
local  field  effects  [47]  are  automatically  taken  into  account 
in  our  calculations  Levine  [32],  [38]  gives  an  expression  of 
the  bond  charge  magnitude  which  g ves  the  best  fit  with 
the  observed  second-order  susceptibilities  of  binarv  com- 
pounds. He  uses  a scalar  model  where  the  different  com- 
ponents a and  at  are  not  Distinguished  For  simplicity, 
we  use  his  result  modified  to  include  the  factor  3 arising 
from  our  tensor  treatment. 

To  calculate  the  nonlinear  suscepnbilitv  of  the  crystal 
we  have  to  add  the  contributions  of  the  different  bonds 


- 97  - 


Cl  IMl«  fl  ,)t  UmD-HtMOMCCIMIrllOS  IN  THI  CMAieOPYR  TI  COO’  At, 

(361.  (41  ] 


i 


75 
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' V A ’•*  A *•» 

Cuki  ■ L c«  0t.  cos 


cos  #*,  ' cos  9,,’  >.,,1  ' (13) 


where  the  index  (s)  covers  the  different  bonds  in  the  unit 
cell  whose  volume  is  V The  capital  indices  U,  J.K.L)  refer 
to  the  crystal  axis(.VKZ),  whereas  the  miner  indices  (i.y.  k. 
/)'  refer  to  the  bond  axis  If.  using  the  axial  symmetry  of 
the  bond,  it  ;s  always  possible  to  choose  the  microscopic 
axis  x ' in  the  X)  pljne.  then  the  angles  appearing  in  (13) 
arc  simple  functions  of  the  polar  angle  \0  * . * * | of  the 
bond  axis  : * . as  shown  in  Fig  2.  The  practical  calculation 
of  (11)  is  then  simplified  by  using  the  operations  of  the 
point  group  of  the  crystal  which  goes  the  directions  of 
equivalent  bonds  in  the  unit  cell 

The  orientations  of  the  two  tvpes  of  bonds  of  a 
clialcopy  rite  compound  .45C,  arc  simply  related  to  the 
geometrical  factors  (36).  (4 1 1 r » 2 - c a and  o * Ax  - I 
If  wC  x,  the  notation  + for  4C  and  - for  BC  to  dis- 
tinguish between  the  two  types  of  bonds,  then  the  relevant 
combinations  of  components  of  r,„  are 

Cm  - 3c,,  * — ™ |(l  — 4(r  + »)  3b 

-F  (I  + 2(r  + „)  9bx*  + (I  - 4(r  - c)  3b,* 

+ (’  * 2(r  - c)  9bi.” I (14a) 


Fif  2 Orxni  jlion  of  the  bond  ant  in  Ihe  crystallofraphic  am. 


TABLE  II 

Comparison  or  tmf  CxrtHivi vial  Third-Order  SusctmaiLiriEt 
Dor  to  tnc  Hound  Fit.  thins  in  Somi  TrTR»i»iDRAi  Semi- 
conductor1; to  Tin  Thioritical  Valuit  Calculated  from 
Our  Mooel 


Compound  r.„  Fxiicrimcniat  c,«  Theoretical 


Gcd..Jc..-(I0±5)I0  "ESU 
1 c .«<6±3)  10  " ESU 


c,, ”S(X  10*"  ESU 
C..-3XI0'"  ESU 


Sj,»«.|c,.-('1«±0  3)  10  " ESUr.,-0? 
c..-(0.3±0  15)10  " ESUc„-0 


9X10  " ESU 


GaAs1 


65X10  " ESU 
[c„-(0  97db0  2)  10  " ESU  c,,»0  95X10  " ESU 
U«-(0  5±0.15)  10  " ESUc„-0.7XI0*"  ESU 


Ci.  ~ c„  - — ((r  + <rb4*  + (r  - cb.'l  (14b) 

r*«  = ! I * ~ (r  -f  #)/3h<*  + 

+ (I  — (r  — <r).  3h  + 2y/ I . (14c) 

Since  the  chalcopyrit;  structure  is  a superstructure  of 
diamond  and  sphalerite  structures.  (14)  can  also  be  used 
for  elementary  and  b.nary  compounds  with  these  struc- 
tures by  putting  r = c * 0. 

The  expression  of  the  deviator  (36)  component  (c„  - 
<■„)  is  to  be  compared  to  another  dev  lator  (36)  fx  - x ). 

• c . the  difference  between  the  parallel  and  transverse 
component  of  the  first-order  susceptihilitv  As  expected, 
the  two  dcviators  exhibit  exactly  the  same  geometrical 
dependence.  Since  r and  <r  arc  small  (10  ' to  10  ')  the 
combination  (r„  — c„)  is  one  or  two  orders  of  magnitude 
smaller  than  the  individual  components  whatever  the 
values  of  y * and  y are 

In  Fable  II  the  calculated  third-order  susceptibilities 
due  to  the  bound  electrons  arc  compared  to  the  ex- 
periment <1  determinations  for  some  'enticandactorv  m 
which  this  contribution  has  been  i>oijted  (2S).  |'(*j.  and 
the  agreement  is  good.  In  the  case  of  CdGcVs,  the  model 


predicts r,,  - 1.8  x 10*"  ESU:c„  - I.,  X IO*l:ESU;and 
- 1.4  X I0-"  ESU. 

Our  analysis  also  pr  ivides  important  information  about 
the  physical  origins  of  the  nonlinear  effects  in  semiconduc- 
tors. An  alternative  of  the  bond-charge  model  (32)  is  the 
charge-transfer  model  (40).  where  the  variation  of  the 
linear  susceptibility  is  attributed  to  a transfer  of  electric 
charge  from  one  atom  to  another  under  the  applied  field. 
The  two  models  describe  the  second-order  susceptibility 
with  comparable  accuracy  But  several  difficulties  arise  in 
the  description  of  the  third-order  susceptibility  when  use  is 
made  of  the  charge  transfer  model. 

First,  there  is  no  transfer  of  electric  charges  when  the 
field  is  perpendicular  to  the  bond  ans  and  y . vanishes  As 
a consequence  c,,  * c„  for  all  the  sphalerite  compounds  in 
disagreement  with  the  experimental  results  (28).  (30) 
Furthermore,  as  only  the  heteronolar  part  of  the  mean 
energy  gap  C is  charge  dependent,  there  is  no  covalent  or 
mixed  contributions  to  the  nonlinear  susceptihilitv  Accord- 
ing to  (12)  y is  proportional  to  (4/;  - I)  leading  to  very 
small  nonlinear  susceptibilities  for  crvstal  with  an  lomcitv 
/,  *■  0 25  For  example,  in  the  ease  of  gallium  arsenide /,  * 
0.27.  the  modcl-prcdict*  : |0  "ESU 
Th:  diserepjncv  'huui ! rot  he  .ittnhuted  to  the  charge* 
transfer  model  itsclt  since  it  :\  well  known  tnat  deforma- 
tion of  the  electronic  distribution  due  to  a field  perpen- 
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dieular  to  the  bond  axis  exists  But  unf  rtunatclv.  the 
formalism  of  Phillips  theors  is  unable  tc  describe  those 
effects  as  well  as  the  contribution’s  second  order  in  7« 
and  7.  when  the  field  is  parallel  to  the  bond  axis 

III  Theory  of  the  Third-Order  SLSCErTtatLm 

or  the  Free  Carriers  in  a Chalcopyrite  Crystal 

The  other  system  of  charges  that  has  been  recognized  to 
make  a large  contribution  to  the  third-order  susceptibility 
of  semiconductors  is  the  free-carrier  gas  There  are  two 
distinct  mechanisms  which  arc  responsible  for  their  non- 
linear susceptibility  the  nonparabolicits  of  the  energv 
bands  Jl3),  (21).  [22]  (i.e  an  energy-dependent  effective 
mass),  and  an  energy-dependent  relaxation  time  (23|-[2R| 
The  two  mechanisms  have  a vcrv  different  frequency 
dependence  The  second  is  only  important  w hen  tlir  inverse 
of  the  relaxation  time  is  not  negligible  compared  to  the  fre- 
quencies of  the  interact  mg  fields  This  can  occur  in  the  case  of 
mixing  experiments  [27]-[29|  where  w,  * u,  -*•  u,  - u,  and 
the  difference  of  frequencies  w,  - w,  is  very  small  flO'*s-'). 
When  only  addition  of  frequencies  is  considered  as  m THG 
(3m  * w + m + w)  the  effects  of  carrier  scattering  become 
negligible  (26]  and  onlv  the  contribution  of  the  non- 
parabolieity  of  the  energv  bands  need  be  considered 

For  elementary  and  binary  compounds,  the  Kane  model 
(31)  gives  a good  description  of  the  energv  bands  and  has 
been  sucessfully  used  to  explain  the  third-order  suscep- 
tibility of  some  semiconductors  (13).  (2l]-[26).  Kildal  has 
recently  extended  the  Kane  model  to  ternary  compounds 
with  ehalcopvnte  structure  (31).  The  principal  feature  of 
Kildal's  extension  is  that  in  addition  to  the  spin-orbit  split- 
ting A one  hcs  to  consider  the  crystal-field  splitting  6 due 
to  the  tetragonal  compression  of  chalcopy  rite  compounds 
A*  a result  of  the  two  perturbations,  the  band  structure  of 
the  ehalcopvnte  ll-IV-V,  crystal  (i  e..  containing  no  noble 
metal)  consists  of  one  conduction  band  (j-like)  and  three 
split-valence  bands  (/»-like)  The  four  bands  ean  be 
described  bv  four  functions  £„’(*)  ■ Emik)  (/  2 m)k'. 
where  the  Em(k)  arc  the  solutions  of  the  fourth-order 
equation 

£(£  - £,  - £eX£  - £,X£  - £,)  - »’,’£(£  + JA) 

- »Fa’[(£  - £,X£  - £,)  - j (£  + i)J  - 0 (15) 


COwAO-AO  CX*.C0»-»|TC  C0w»3uh? 


ro«  CiC4t*i  [(  ODiv  0 ■ OJliv  l-  oniv 
t,  -OMiv  (,  -OXl* 

Ft?  J Schematic  of  the  hind  oructure  of  a chitcnpvrile  compound  xod 
iu  ipholeriie  analog 


transitions,  the  third-order  susceptibility  due  to  the  non- 
parabolicitv  of  the  energy  band  is  given  by  [13],  (21).  |22] 


i - 3c,....  ; : — 

*4  n T*  + W,) 

. y i if  

4*  ,r  ik.Akfk^k,  (16) 

where  /»(£)  is  the  Fcrmi-Dirac  distribution  function 
For  an  n-type  material  the  positive  sign  and  Er  arc  used, 
whereas  for  a p-tvpe  compound  the  negative  sign  and  £v , 
arc  used 

For  the  elemcn  ary  and  binarv  compounds  when  the 
energy  gap  £r  is  small  enough  so  that  the  deep  valence 
band  £>,  is  unperturbed,  a closed  form  can  be  deduced 
from  Kane’s  model  (33)  which  gives  a good  description  of 
£r  and  £v , This  close  form  has  been  widely  used  for  com- 
pounds with  very  small  band  gap  like  InSb  For  Crystals 
*vuh  a larger  band  gap  the  closed  form  is  not  valid,  and  a 
better  description  af  the  band  shapes  i«  given  by  expand- 
ing the  energv  as  a funr'ion  of  the  wave  vector  [43]  In 
the  same  line  we  assume  the  following  form  for  the  energy 
bands  of  a ehalcopvnte  compound 

£-(*)-£..  ±AmW  * ♦ cmH  *+  D„H\* 

+ Gm  H’»”+  (17) 


Here  W * - P>k,\  - £*(*,«  + A ,*)./»-  - < 

S[p,\Z  > [29],  (30],  and  £,  and  £,  arc  solutions  of  £,*  + 
(a  + A)  £|  + (2Aa/3)  - 0 

In  the  limit  b ■ 0.  ( 1 5)  reduces  to  the  Kane  equation  The 
band  structures  of  a chaldopyritc  compound  and  of  its 
sphalerite  analog  arc  shown  in  Fig  3.  Fig  3 also  lists  the 
values  of  the  parameters  deduced  bv  Kildal  [31]  to 
describe  the  band  structure  of  CdGcA« 

W’hen  the  frequencies  'f  the  inieracii  field  are  small 
compared  to  the  frequencies  associated  with  the  electronic 


In  this  equation  the  index  m runs  over  the  four  enerev 
band?  The  coefficient',  of  the  expansion  arc  found  bv  in- 
serting (17)  into  (15)  The  results  are  given  in  Table  III 
It  should  be  noted  that  while  the  limits  of  the  Am  for  b = 
0 are  correct,  this  is  not  the  case  of  the  Bm  To  calculate  the 
Bm  one  has  to  divide  bv  £,£,  = (2  *)Ai.  and  this  process  is 
not  valid  for  b or  A caing  to  zero  Therefore,  we  will 
assume  that  in  the  limit  ) « 0 4„  - /?„ 

The  expressions  for  the  coefficients  are  simple,  and  due  to 
the  fact  tl.at  (15)  is  of  the  second  order  in  M and  It  , the 
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TABLE  Ml 

Cot ITIf  It  N 1 1 Of  THt  EXPANSION  Of  Tin  Eniki.v  UaSU  At  A Kmtion 


TABLt  IV 

Pma-j-Maicming  A nclls  ton  THO  and  SHG 


Of  THt  * 

AM  VlClOt 



. 



— 

A, 

Phase- Matching 
C on  figuration 

SHG  1 

SHG  II 

THG  1 

THG  II 

(£,  - £,X£.  - £,) 
At  - 0 

2(£,  - £,X£.  - £,) 

a,- 

*m  (calculated) 
(.  (measured) 

35*  ± 7* 
32*  * I* 

55*  ± 7* 
49*  ± I* 

50*  ± 7* 
(50*) 

67*  ± 7* 
60*  ± 1* 

Ai  ■ 


A.  - 


£,  + U 


(£.  - £,X£,  - £,) 


£«  ± ja 


. _ c»  + 1 + 

J " 2(£,  - £,X£,  - £») 


(£,  - £.X£»  - £,) 


B. 


2 E, 


E,  + I + \* 
HE,  - E.KEi  ~ E. j 


C.  - AmS(Am)  0.  - BmS(B.) 

C.  - AmS(3.)  + fl.SM.) 

where 


SM.) 


E j~zrj~  **.)  - E 


AW  £_(*)  - £.  + A.Vt‘  + *.HV  + C.#V  + fl.*V  + 


The  nonlincarities  due  to  the  free  carriers  are  about  one 
order  of  magnitude  larger  than  the  nonlincarities  due  to 
the  bound  electrons. 


coefficients  satisfy  the  interesting  relation 


IV.  Experiments  and  Discussion 


Am*Dm  + BmlCm  m AmBmGm-  (18) 

Since  we  have  only  retained  the  terms  in  £„  up  to  the 
fourth  power  of  the  wave  vector  the  fourth  derivatives 
appearing  in  (16)  are  constants  and  the  nonlinear  suscep- 
tibilities arcproportionaltothe  number  cfthefreecarrieriV 
The  relevant  combinations  of  coefficients  appearing  in 
the  effective  nonlinear  susceptibility  are 


C|« 


**!•  **!• 


±l8t'(Aw)4 


G. 


V 


5: 


181  (Aw)4 


(G. 


V 


C11  3ci* 

20,  y) 


0 


(19) 


where /«  is  the  energy  of  the  laser  photons 
For  CdGeAs,  one  finds  by  using  Table  III  ard  the 
energies  of  Fig  3 and  £*  - (hV2m)  X 20  eV 


f„*  - -2  5 X I0",,N  r./  - r«*  - 2 X I<T,tN 

c,»  - -1.5  X 10  f,N  c„*  - c„‘  - 6.3  X lO'^N. 

(20) 

The  particular  sample  we  have  used  for  type-ll  PM-THG 
is  a p-type  crystal  with  ,V  * 5 x 10'*  holes/ cm*  so  that 


c„  * - 12  5 x 10  " fcSU  c„  -c„  - 10  x 10  " ESU 


In  this  section  we  present  experimental  work  on  the 
measurement  of  the  third-order  nonlinear  susceptibility. 
CdGeAs,  has  a birefringence  large  enough  (II).  (31)  to 
allow  type-I  and  tvpe-11  PM-THG  for  a fundamental  wave 
length  of  10.6  p.  The  phase-matching  angles  calculated 
(31)  from  the  measured  refractive  indices  (31)  are  listed  in 
Table  IV  for  both  THG  and  SHG  The  difference 
^"(THG)  -t„"(SHG)  * 12*  is  small  so  that  both  PM- 
SHG  and  PM-THG  can  be  observed  in  the  same  sample. 
Thus  the  ratio  of  the  third-order  and  second-order  non- 
linear coefficients  c,,  d„  can  be  measured  with  the  same 
crystal  sample  Furthermore,  at  high  power  densities  a 
substantial  nonphase-matched  mixing  of  the  laser  radia- 
tion with  its  PM-SHG  can  be  observed.  The  comparison 
of  the  two  signa's  at  3w  givrs  a simple  way  to  eliminate  the 
detector  wavelength  response  and  gives  an  accurate 
measurement  of  the  ratio  Cu/d ,«• 

3Ve  have  examined  two  samples  Sr  and  S,,.  Si  was  cut 
for  type- 1 THG  and  5n  for  type  II  The  CdGeAs,  crystals 
were  grow  n at  tne  Center  for  Materials  Research  at  Stan- 
ford University  (49)  The  experimental  setup  is  shown  in 
Fig.  4 The  electrode  structure  of  the  TTA-CO,  laser  is 
very  similar  to  that  described  in  (50)  The  laser  operates  at 
a peak  power  of  50  kW  with  a pulsewidth  of  200  ns  in  the 
TEM*,  mode,  jt  a repetition  rate  between  5 and  25  pulses 
per  second  tpps).  The  laser  peak  output  power  is  stable  to 
better  than  40  percent  over  periods  of  about  1 h An  inter- 
nal Brewster  window  can  be  rotated  to  adjust  the  polariza 
tion.  During  the  measurements,  the  orientation  of  the 
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la^cr  polarization  and  the  generated  harmonics  was 
carefully  cheeked  with  jn  infrared  reflection  polarizer 
In  the  first  set  of  experiments  S,,  was  examined  The 
measured  phase-matching  angles  for  tvpe-ll  TUG  and 
SHG  are  given  in  Table  IV  In  light  of  the  accuracv  ol  the 
refractive  indices  the  agreement  with  the  predicted  values 
is  satisfactory  It  should  be  noted  that  there  is  a better 
agreement  on  the  difference  0„"(THG)  - IJ'fSHG) 
The  signal  at  3w  from  the  direct  THG  and  from  thesum 
generation  have  a very  different  behavior  The  angular 
phase-matching  tolerance  is  very  acute  on  the  THG~  and 
the  half  width  of  the  phasc-m  itching  peak  .50  = 10*  t I* 
(outside  the  crystal)  is  in  good  agreement  with  the 
predicted  value  SO  - 12*  ± 2*  for  the  I 75-mm  crsstal  As 
the  crystal  is  rotated  the  signal  at  3w  goes  to  zero  then  in- 
creases and  has  a very  broad  maximum  around 
0„"(SHG).  which  can  be  maximized  bv  translating  the 
slightly  wedged  sample  By  selective  use  of  a sapphire  and 
fuscd-silica  filter  we  measure  the  SH  and  TH  signal  with  a 
liquid-nitrogen-cooled  InSb  detector 
The  laser  is  focused  with  a germanium  lens  such  that  the 
aperture  length  (5 1 ) /,  ■ 17  mm  and  the  length  of  the  focus 
(52)  I,  » 220  mm  are  much  larger  than  the  sample 
thickness  / ■ I 75  mm  The  plane-wave  analssis  is 
therefore  valid  (53),  and  the  harmonic  power  at  phase 
matching  can  be  calculated  according  to  (54).  To  calculate 
the  intensity  at  3w  due  to  sum  generation  we  assumed  that 
the  field  at  2 w-  is  constant  over  the  last  coherence  lencth  /r 
and  has  the  value  - ( I '2)1  ).  and  we  assumed  that  the 
field  at  3w  is  mostly  radiated  in  a thin  slab  of  material  with 
a thickness  lt  at  the  output  surface  of  the  crystal  The 
powers  at  2u  and  3w  anj  2w  w are  then  given  bv 
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where  w,  and  n,  are  the  beam  waists  and  the  absorption 
coefficients  at  the  frequenev  nw.  t is  the  transmission  fac- 
tor for  the  laser  field  at  the  input  surface  of  the  crsstal.  r, 
and  r,  the  transmission  factors  for  the  harmonic  fields  ,u 
the  output  surface.  and  c,,.  d,,.  and  du  arc  the  effective 
nonlinear  coefficients  for  THG.  SHG.  and  nnxme  The 
angular  dependence  of  the  transmission  factors  has  to  he 
taken  into  account  'ince  the  incidence  angle  lor  PM-SIIG 
is  quite  large  (>43*).  The  absorption  coefficients 
measured  for  the  sample  with  a CO,  laser,  its  SH.  and  with 
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an  He-Ne  laser  at  3 '9  „ are n,  * 0.5  cm  '*.«,»  1 .3  cm  !. 
<»»  * 5 7 cm  '. 

The  measurement  of  THG  to  the  sum  generation  leads 
to  IrJ  * (170  i 251c/,.1  To  compare  the  THG  io  the 
SHG  we  have  to  consider  the  frequencs  dependence  of  the 
sensitivity  of  the  InSb  detector  We  estimate  o’,r  particular 
detector  to  be  twice  as  sensitive  at  5 1 u as  at  3 s ^ This 
gives  a measured  vahe  of  IrJ  - (|  6 ± 0 8)10  V«!  bv 
comparison  of  PM-ThG  to  PM-SHG  If  wC  use  the  value 
85  x 10  ’ ESU  for  th-  d„  (see  Appendix  II  for  dis- 
cussion), the  two  measurement  *n  good  agreement 
and  give,  respectively.  Ir,,|  * 12  2 x 10  " LSU  and  If,,) 

13  6 X 10  " ESU  In  light  of  the  accuracv  of  the 
measurement  and  of  the  value  of  «/„  we  obtain  |r„|  * ( 1 1 
± 6)10  " ESU. 

In  the  second  set  of  experiments  we  compared  PH-THG 
in  Si  and  Sh  The  experin  ents  were  difficult  to  perform 
cue  to  the  much  smaller  magnitude  of  the  type-1  PM-THG 
signal  Also  mixing  of  the  laser  with  the  nmi-PM  SHG 
produces  a faint  signal  at  3 5 p which  could  be  observed  at 
all  angles  of  incidence  Furthermore,  the  phase  matching 
for  type-1  THG  and  for  type-ll  SHG  are  ver;.  close,  and  a 
permanent  use  of  the  IR  polarizer  was  necessary.  The 
recorded  signal  at  3«e  showed  a broad  maximum  around  0 
* 50°.  but  the  signal-to-r.oise  rano  was  too  poor  to  surelv 
assign  our  observation  o type-1  THG  alone  However,  an 
upper  limit  to  the  tvpe-l  THG  coefficient  was  deter- 
mined to  be  If, i ^ 5 «■  10  J lr,, | * 06  * 10  " ESI 
Using  the  theoretical  results  of  Sections  II  and  III  and 
(3)  for  the  expression  ol  the  effective  nonlinear  coefficients 
gives,  for  the  theorct  ;al  values,  c,  * -0.15  X 10  :l  ESU 
andfi,  * — 4x  10  ESU.  The  tvpe-ll  effective  nonlinear 
coefficient  which  i«  mainlv  due  to  the  free-carrier  contribu- 
tion is  between  two  and  four  times  larger  than  predicted 
from  the  Kane  model  Such  a difference  was  a Ireadv  noted 
for  n-tvpe  GaAs  bv  Cardona  )4SJ  m measurements  of  the 
reflectivity  and  .he  I arudav  -otatior.  as  a function  of  the 
number  of  carriers,  and  bv  Wvnnc  [29)  in  measurements 
of  the  third-order  susceptibility  in  mixing  experiments 
These  authors  also  fou,id  that  the  nonparabolicitv 
described  bv  the  cosvF  icnts  ol  x*  in  the  expansion  of  the 
energy  was  larger  than  predicted  hv  the  Kane  model,  (two 
times  for  the  results  of  Cardona  and  three  times  for  those 
of  Wynnel  It  should  be  noted  that  since  the  effective 
masses  of  the  holes  on  6,  are  quite  large  (3 1 ) w * 0031  n 
and  m - 0 m.  the  For*-”  energy  in  our  sample  is  sma1!. 
£p  * 2 meV.  compared  to  the  splitting.  200 

meV  Therefore,  no  additional  contributions  Irom  the 
lower  valence  hands  /,.  and  £,,  can  explain  the 
observed  difference 

Reversing  the  above  areument.  we  can  sav  that  the 
measured  value  of . ,,  is  conostent  with  a r.onparabolicitv 
of  the  conduction  bard  three  times  larger  than  predicted 
from  the  Kane  model  \s  tor  . , the  theoretical  value  s 
al<o  smaller  than  the  measured  one  The  sigrjl  that  we 
observed  should  theretorc  be  .'ltrihutcd  to  sain  ecncr.i- 
tion  The  sample  used  was  verv  thin  tl  mint  ,nd  when 
larger  erys*,t|>  .,re  available  i 's  4 nnn)  it  should  be  possib  e 
to  observe  ivp»-l  PM  TUG 
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C Ml  Ml  « rl  ul  TMI*S)  H;,»MOSir  Of  Ntll.THlN  in  rill  cm.i  1 l>n  mil  MV.t.S 


79 


F inally,  we  would  like  to  discuss  Sriefls  the  practical  use 
of  CdGcA  Sj  as  an  IR  harmonic  generator  Due  *o  the 
large  value  of  the  coefficient  ct|  and  the  small  v alue  of  the 
difference  0„"(THG)  - 0„"(SHG),  high  conversion 
efficiency  can  be  obtained  for  both  SHG  and  THG  in  the 
same  sample  For  example,  using  a mode-locked  CO,- 
Tr  A laser  ofO  I J/pulsc.  a pulse  duration  of  I ns.  and  a I- 
cm-long  crystal  it  ts  possible  to  reach  30-perccnt-THG  and 
60-perecnt-SHG  conversion  efficiency. 

IV.  Conclusion 

We  have  analyzed  the  possible  THG  configurations  in 
ch  tlcopyritc  crystal  from  the  point  of  view  of  the  irreduci- 
ble components.  The  bound-electron  contribution  and  the 
free-carrier  contribution  to  the  third-order  susceptibility 
can  be  expci  mentally  separated  bv  the  appropriate  THG 
experiments  We  have  generalized  the  bond-charge  model 
to  calculate  the  nonlinear  susccptibilitv  of  the  bound  elec- 
trons The  results  are  in  good  agreement  with  the  ex- 
perimental determinations.  We  have  used  Kildal's  exten- 
sion of  the  Kane  model  to  chalcopyritc  compounds  to 
calculate  the  free-carrier  contribution  ° *r  measurement 
on  CdGcAs,  shows  that  the  nonparabolicitv  of  the  first 
conduction  band  is  about  three  times  larger  than 
predicted.  The  third-order  nonlinear  coefficient  of 
CdGcAs,  has  a large  value  c,i  - 12.5  y 10’ " ESU  so  that 
very  high  conversion  efficiency  for  both  SHG  and  THG 
can  be  achieved  i i this  crystal. 

Atpendix  I 

Contracted  Notation  for  THG 


ARPENDtX  II 

Second-Order  Nonlinear  Coefficient  of  CdGeAs, 

1 he  second-order  nonlinear  coefficient  of  CdGeAs,  has 
been  measured  relative  to  J.,  (GaAs) 

According  to  Byer  rial.  |3). </„(CdGeAs,)/<MGaAs)  - 
3 4 ± 25  percent.  From  Boyd  et  al  (4) 

«/i«(CdGeAs;)/</M(GaAs)  * 2 6 ± 15  percent  We  use  the 
value  <UCdGcAs,)  * 3d„(GaAs) 

The  </, 4 coefficient  of  GaAs  has  been  measured  com- 
paratively to  Ag,SbS,  and  CdSe  whose  coefficients  have 
been  measured  absolutely 

</,.(GaAs>  - (7  5 ± 0 3)  i/.(A*,SbS,)  |55| 
t/ufCdSc)  - (0  2 -t  0 02)  tUGaAs)  (2] 

and 

<!  fAg,SbS,)  ■»  0 43  10  ESU  ± 30  percent  |56J 

f/,,(CdSe)  ■»  O 53  10  ESU  ± 15  percent  |57) 

These  sets  of  values  fit  well  with  </,.(GaAs)  - 2.8  X '0  ’ 
ESC  ± 40  percent,  and  therefore  give </„(CdGeAs,)  * 8 5 
x 10-’  ESU. 
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Throughout  this  article,  we  used  the  notation  of 
Midwinter  and  Warner  [20J.  The  contracted  notation  is 
defined  by 

w I 2 3 4 5 6 7 8 9 10 

ijk  III  222  333  233  223  133  tl3  122  112  123. 

The  nonlinear  polarization  at  the  TH  frequency  is  given  by 
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Second  harmonic  generation  and  infrared  mixing  in  AgGaSe, 

R.  L Byer,  M.  M.  Choy,  R.  L Herbst,  D.  S.  Chemla*,  and  R.  S.  Feigelson' 

Microwatt  Laboratory.  W W Haiuen  Laboratories  of  Physics.  Stanford  Unirersity.  Stanford,  California  94)05 
(Received  17  September  1973;  in  inal  form  29  October  1973) 

We  have  continuously  tuned  between  7 and  I)  p.m  by  mixrt|  the  output  of  a LiNbO.  parametric 
oscillator  tn  the  chalcopyrtte  AfOaSe  We  have  doubled  a CO,  laser  with  2 7%  efficiency  which 
agrees  very  well  with  the  expected  efficiency  and  verifies  the  high  optical  quality  of  the  I 'J-cm-lon* 

AgGaSe,  crystal.  The  measured  transparency  range,  indices  of  refraction,  and  nonlinear  coefficient  of 
d n m (3  7 £ 0 6)  x 10  11  m/V  show  that  AgGaSe,  is  a useful  infrared  nonlinear  material  phase 
matchable  over  the  entire  J-18-pm  infrared  region. 


Since  the  first  demonstration  of  phase -matched 
second  harmonic  generation  (SHG)  in  AgGaS,, 1 the  non- 
iinear  properties  of  the  ternary  semiconductors  with 
chaicopyrite  structure  have  been  wideiy  studied.*”* 
Their  iarge  nonlinear  susceptibilities  together  with  ade- 
quate hirefringence  to  achieve  phase  matching  make 
them  attractive  for  nonlinear  optical  devices.  Nonlinear 
mixing  has  been  demonstrated  in  ZnGePb  ’ AgGaS2,  *'* 
CdGcAs*, l0,u  and  recently  AgGaSe,. M,u 

AgGaSe,  singie  crystals  are  grown  by  the  vertical 
Bridgeman  method  after  the  starting  materials  are 
presynthesized  in  a carbon  boat  contained  in  a seaied 
quartz  crucibie  The  presynthesized  stoichiometric  mix 
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with  a meiting  point  i approximately  660  C is  then 
transferred  to  a quartz  crucible  coated  with  pyroiiGc 
carbon  for  a pregrovth  run  at  2 mra/h  rat'1  through  9 
40  C/cm  temperature  gradient  The  top  and  «.Mto’a  of 
the  resuiting  bouie  are  then  removed  prior  to  the  actual 
growth  run  in  the  same  vertical  furnace  at  a siow  growth 
rate  of  0. 2 mm  A.  After  growth  the  crucible  is  cooied 
to  room  temperature  at  25  “C  A-  The  resuiting  14-mm- 
diam  single  crystals  typically  show  a galiium-rich  re- 
gion near  the  seed  end  foiiowed  by  an  approximately  2 
cm  useful  AgGaSe,  singie -crystal  region  and  a siiver- 
ruch  top  section.  Eariy  crystals  6howed  a precipitate 
which  reEuited  in  a 2-cm*'  scatter  loss.  In  recent  crys- 
tals the  scatter  ioss  has  been  significantly  reduced  to 
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Second  harmonic  feneration  and  infrared  rmxmf 
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FIG.  1.  .Tie  mixing  output  <Xlr>  in  angle  phase- 

matched  AjGaSe,  pumped  by  a 1.318-pm  Nd:  VAG  liner  mix- 
ing with  the  output  of  a 0.  o59-pm  pumped  LINbO,  parametric 
oeclllntor  (Xlx). 


near  0 02  cm*1  at  10.6  pm.  All  crystal!  have  high 
resistivity  and  good  optical  transparency  from  the  band 
Kap  at  0. 71  pm  to  the  taro-phonon  absorption  edge  at  16 
pm.  A particularly  attractive  feature  of  AgGaSe,  is  the 
case  with  which  single  crystals  can  be  grown. 

Wc  have  measured  the  nonlinear  coefficient  of 
AgGaSe,  relative  to  GaAs  at  10. 6 pm  using  the  wedge 
technique  u~u  The  wedge  method  is  suited  for  simill- 
blrcfringence  large -index  materials.  For  the  wedge 
method  the  coherence  length  is  given  by 

~nM)  = ,Ay  tana, 

where  Av  is  the  wedge  translation  distance  between 
SHG  extrema,  a is  the  wedge  angle,  and  the  beam  is 
'ncldcnt  normal  to  the  input  face  of  the  wedge.  For  our 
crystal  samples  the  coherence  lengths  are  expected  to 
bo  near  100  pm  for  GaAs  and  200  pm  for  AgGaSe,  at 
10.6  pm. 

Wo  used  a (OQl)-oriented  AgGaSe,  wedge  for  the  non- 
linear susceptibility  measurement  relative  to  two  GaAs 
reference  samples.  One  (001) -oriented  wedge  was  cut 
from  GaAs  material  supplied  by  Coherent  Radiation 
Laboratories  and  a second  (111) -oriented  wedge  was  cut 
from  Monsanto  material;  both  GaAs  samples  were 
chromium  compensated  to  a high  resistivity  of  about 
10  Ocm.  With  the  laser  polarization  parallel  to  the 
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(110)  direction,  £,=  £,=  £*//? and  P=P,=  €*#„£*,  so 
the  rf„,  for  the  (001) -cut  GaAs  and  AgGaSe,  samples 
are  du  and  dM  = dx„  respectively.  For  the  (lU)-cut 
GaAs,  with  the  laser  polarization  along  tbs  (111)  direc- 
Uon-  E>=  E.  = EJ^ and  Pt  = Pf  = pt  = p/4 \ so 
2/V3 d14Eg  and  the  d«t  In  this  case  is  2//3 ~du  The 
coherence  length  of  GaAs  is  measured  to  be  107*  2 pm 
and  109* 3 pm  for  the  (001)-  aid  (lll)-cut  samples. 
This  is  in  excellent  agreement  with  published  values  In 
Refs.  17  and  18  of  107*5  pm  and  107*  1 pm.  The  mea- 
sured coherence  length  of  AgGaSe,  Is  237*  15  pm  which 
compares  favorably  with  a calculated  value  of  255*  50 
pm  based  on  the  index  data  of  Boyd  et  at. ' 

The  AgGaSe,  nonlinear  coefficient  measured  relative 
to  GaAs  is 


(AgGaSe,  )/du(GaAs)001  = 0. 33  * 25%, 

= rf*(AgC.aSe,)/rfu(GaAs)  1 1 1 = 0. 32  * 18%. 

These  values  are  in  good  agreement  with  the  relative 
measurement  by  Boyd  et  at.*  of  rfM(AgGaSe,)/rfM(GaAa) 
*0.37*10%.  Taking  du(GaAs)=  (11. 7*4C%)xio-»  m/ 
V,  we  obtain  dM(AgGaSe,)=  (3  8*  1. 7)xlorn  m/V. 


We  have  also  performed  phase -matched  SHG  of  10  6 
urn  using  a 60  -cut  AgC-Ve,  crystal.  The  measured 
phase -matching  angle  of  57.5°*  0.5*  is  In  good  agree- 
ment with  the  calculated  value  of  55*  *4*.  The  expected 
phase  matched  SHG  conversion  efficiency  is 

Pu/P»=  rV=  (2u)ld£, /wnin^ (B,  f),  (2) 

where  the  powers  are  defined  inside  the  crystal 
- d sin0^,,  / is  the  crystal  length,  k^2m^/x,  and 
h\B,  f)  is  the  Boyd  and  KleinmanM  focusing  factor  which 
reduces  to  l/b  = l/u>lk  in  the  loose  focusing  limit.  For  a 
low-loss  AgGaSe,  crystal  1 cm  in  length  In  the  loose 
focusing  limit  (l/b<  1)  the  calculated  SHG  conversion 
efficiency  is  r */*=  0. 75%  at  1 MW /cm'.  Using  a TEA 
CO,  laser  operating  in  a TEM*,  mode  as  a source,  we 
mearired  the  absolute  SHG  efficiency  generated  in  a 
high-quality  1. 54-cm-long  AgGaSe,  crystal.  The  aver- 
age input  and  output  powers  measured  with  an  Eppiey 
thermopile  were  2. 82  mW  and  76  pW,  which  corre- 
sponds to  1. 6 kW  of  peak  power  at  the  fundamental  for 
a 200-nsec  pulse  width.  The  experimentally  observed 
conversion  efficiency  for  the  incident  Intensity  of  1. 68 
MW/cm  is  2.63%  Using  38 xl0"u  m/V,  t 54 
cn  , and  h(B,  £)  = 0 925  l/b,  which  corresponds  to  the 
focal  spot  size  of  250  pm,  the  expected  conversion  ef- 
ficiency is  2.76%.  This  measurement  can  be  considered 
*s  a separate  absolute  determination  of  the  nonlinear 
coefficient  of  AgGaSe,.  The  nonlinear  coefficient  is 
found  to  be 


^h=  (3. 74  * 0.  6)  >.l(Tn  m/V, 

which  agrees  very  well  with  the  previous  measurement 
made  relative  to  GaAs  and  with  the  recent  absclute 
measurement  of  Kildal  and  Mikkelsen"  of  dM=  (3. 24 
*0.  50)xl(Tn  m/V. 

AgGaCc,  prese  matches  for  SHG  for  fundamental 
wavelengths  between  3 and  13  pm.  * The  SHG  efficiency 
of  AgGaSe,  is  significantly  better  than  proustite,  for 
example,  due  to  both  a factor-of-three  increase  in  the 
nonlinear  coefficient  and  the  small  birefringence  which 
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FIG.  2.  Phase-matched  mixing  peak  generated  in  a fixed 
AgGaSe,  cy ratal  pumped  with  a tunable  source. 


allows  increased  interaction  lengths  without  aperture 
length  limitations. 

AgGaSe,  has  adequate  birefringence  to  phase  match 
over  an  extended  infrared  spectral  range  for  tunabie 
wavelength  generation  by  mixing.  For  mixing,  a con- 
venient tunabie  pump  source  is  a LiNbO,  parametric 
oscillator11  mixing  with  a fixed-frequency  Nd:  YAG 
laser.  For  our  experiment  we  used  a collinear  geome- 
try. The  acoustic  Q-switched  Nd  YAG  laser  tuned  to 
1.32  pm  Is  internally  doubled  with  a LilO,  crystal  to 
generate  0.659  pm  which  pumps  the  temperature -tuned 
LiNbO)  oscillator.  The  oscillator  output  in  the  1.5  — 
1.7-pm  range  mixes  in  the  AgGaSe,  crystal  with  the  re- 
maining collinear  1.32-pm  beam  which  passes  through 
the  parametric  oscillator.  The  AgGaSe,  is  angle  phase 
matched  by  rotation  on  a geared  stage. 

The  mixing  efficiency  is  given  by 

PjP-,=  (wir/w-JrVslnc^AM  /2),  (3) 

where  T *!*  — (2 w,ra’M//*(l/*/n>»#,^(lrr,c0)/>  and  lp  is  the 
pump  intensity.  For  a 1 -cm  AgGaSe,  crystal  and  a pump 
Intensity  of  1 MW /cm*  at  1.32  pm,  the  conversion  ef- 
ficiency is  P|r/PMe=  1.2^(u>,r/u>Mt). 

Figure  1 shows  the  generated  mixed  output  from  7 to 
15  ptr.  Beyond  15  pm  our  HgCdTe  detector  is  response 
limited.  Figure  2 shows  the  phaie  matching  peak 
generated  by  mixing  in  a fixec  AgGaSe,  crystal  pumped 
by  a tunabie  LiNbO,  parametric  oscillator.  The  charac- 
teristic sinc<(AM/2)  phase -mate hing  peak  width  agrees 
with  that  calculated  from  the  dispersion  of  AgGaSe,. 

The  24  6-cm'1  width  is  the  acceptance  bandwidth  of  the 
mixing  crystal  which  is  much  greater  than  the  output 
bandwidth  of  2 cm'1  determined  by  the  2-cm'1  gain  line- 
width  of  the  parametric  oscillator. 

A piot  of  the  phase -matching  rngle  versus  LiNbO, 
oven  temperature  showed  a nearly  linear  relation  over 
a wide  7 — 12-pm  spectral  range.  We  therefore  used  a 


stepping  motor  and  synchronously  rotated  the  AgGaSe, 
crystal  to  phase  match  with  a 1 C/min  temperature 
scanned  LiNbO,  parametric  oscillator.  In  this  wav  the 
spectrum  between  7 and  12  pm  was  continuously  tuned 
8 min. 

Figure  3 shows  a spectrum  of  polystyrene  as  an  ex- 
ample of  the  continuous  scanning  capability  of  this 
source.  The  spectrum  > is  taken  using  a dual  channel 
differential  boxcar  with  *o  HgCdTe  detectors,  Th« 
mixed  c jtput  is  detected  with  better  thin  a 30-dB  signal - 
to-noise  ratio  with  a peak  to  peak  va  iation  of  less  than 
ltfl  at  a repetition  rate  between  10  and  25  pps 

When  mixing  against  1.32  pm,  AgGaSe,  does  not  have 
adequate  birefringence  to  phase  match  at  wavelengths 
shorter  than  7 pm.  Based  on  phase -matching  calcula- 
tions using  the  index  of  refraction  data  of  Boyd  el  at.,* 
mixing  against  w velengths  1.  5 pm  and  longer  allows 
complete  coverage  of  the  infrared  As  an  example,  a 
1.06-pm  pumped  LiNbO,  parametric  oscillator  with 
degeneracy  near  2. 12  pm  angle  tunes  over  a 1. 5-3.7- 
p:n  range.  AgGaSe,  phase  matches  for  mixing  the  signal 
and  idier  waves  to  generate  3-18  pm  for  phase  match- 
ing angles  between  80r  and  50  . u This  example  shows 
the  unique  phase -mate hing  properties  of  AgGaSe,  for 
extended  infrared  generation  by  mixing. 

In  conclusion,  we  have  measured  the  nonlinear  coef- 
ficient of  AgGaSe,  and  demonstrated  phase -matched 
SHG  of  a CO,  laser  as  a verification  of  crystal  quality 


FIG.  3.  Spectrum  of  polystyrene  taken  with  ta)  the  AgGaSe, 
mixer  (resolution  2 cm"1)  and  tbi  Perkin  Elmer 
spect  rophotomete  r. 
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and  potential  use  as  a second  harmonic  generator.  Us- 
ing a LiNhOj  parametric  oscillator  as  source,  we  have 
generated  continuously  tunable  output  between  7 and  15 
pm  by  mixing  in  AgCaSe,. 
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ABSTRACT 

The  Bond  Orbital  Model  for  tetrahedral  compounds  is  used  self- 
consistently  to  calculate  the  second  order  susceptibility.  No  adjustable 
parameters  are  used  and  agreement  with  experiment  is  good.  The  origin  of 
the  nonlinearity  arises  from  a charge  transfer  as  a result  of  the  asym- 
metric polar  energy  between  the  anion  and  the  cation.  The  model  cor- 
rectly predicts  an  optimum  polarity  to  maximize  second  order  suscepti- 
oilities. 


On  leave  of  absence  from  The  Middle  East  Technical 
Turkey. 


University, 
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HJTRODUCTICN 

The  understanding  of  the  origin  of  the  nonlinear  optical  susceptibilities 
Las  been  the  subject  of  zany  'investigations.  This  understanding  is  important 
* in  material  assesnen:  and  prediction  of  new  material  characteristics  in  non- 
linear optical  devices  such  as  harmonic  generators,  mixers  end  parametric 
oscillators.  * In  a more  formal  approach,  quantum  mechanical  expressions 

•re  derived  from  per1  orbation  expansion  involving  complex  sums  of  matrix 

o 5 

elements  and  energy  denominators. * These  fundamental  expressions,  even 

vber*  simplified  through  the  use  of  octupole  moment  of  the  ground  state  charge 

6 7 

distribution,  or  by  employing  ground  state  wave  functions,  are  still  formidable 

to  use  in  quantitative  predictions.  Exact  wave  functions  in  solids  are  often 

not  known  to  high  accuracy  and  the  complete  band  structure  is  also  required. 

Other  complications  from  contributions  to  the  nonlinearity  from  different 

8 9 

points  in  the  Brillouin  zone  can  also  arise.  * Another  approach  is  to  use 

simple  phenomenological  models  such  as  the  enharmonic  oscillator  mcdel^  and 
••  • 11 

Killer's  rule.  Included  in  this  category  are  the  universal  semiconductor 

12  11  14  15 

model,  * the  bond  charge  model,  and  the  charge  transfer  model,  and  the 

16 

molecular  orbital  model.  The  latter  three  models  all  draw  heavily  from 
the  results  of  Phillip's  dielectric  theory. 

‘ 19  2C 

In  this  paper  we  use  an  approach  based  on  Harrison's  3cnd  Orbital  Model.  ’ 

This  is  a more  predictive  model  and  only  appeals  to  essentially  two  input  factors, 

the  dielectric  constant  and  the  E„A(X^-X^)  transition  energy  to  calculate  the 

model  parameters.  The  model  predicts  useful  properties  related  to  valence 

21 

electrons  such  as  cohesive  energy,  band  structure,  etc.  Using  3CM 
parameters,  we  calculate  the  second  order  susceptibility  in  the 


111 


optical  transparency  region  of  a crystal.  Th$  analysis  clearly  displays  the 
physical  origin  of  the  nonlinearity*  Perhaps,  more  importantly,  by  expressing 
the  nonlinearity  in  terns  of  the  fundamental  material  parameters  of  polarity  and 
the  bond  length,  salient  characteristics  desirable  for  high  optical  non- 
linearities  are  distinctly  revealed. 

II.  BOND  ORBITAL  MODEL 

In  this  section,  we  outline  the  BOM  so  that  we  can  discuss  its 
extension  to  the  calculation  of  the  second  order  susceptibility.  For  a 
■ore  complete  discussion  of  the  model  and  the  implications  of  the  approximations 

Involved  refer  to  the  original  papers.  ^*20 

The  model  is  based  on  a tight-binding  approximation  of  sp^  hybrid 

orbitals  for  tetrahedral  compounds.  Denoting  a sp^  hybrid  orbital  on  an  anion 
(cation)  by  |h*  > (|hC>)  , a bonding  orbital  is  constructed  from  a linear 
combination  of  nearest  neighbor  hybrids  on  the  anion  and  its  adjacent  cation 
directed  along  the  sane  bond, 

|\,  > - u#  |h*  > ♦ «c  |hC  > (1)  ' 

The  minimization  of  the  expectation  value  oi  *be  hamiltonian  by  a standard 
variation  on  the  ue  and  ufi  results  in  a (2  x 2)  secular  equation.  The 
solution  of  the  equations  gives  the  eigenvalue  and  the  eigenvectors.  They 

ere 


Vo  < 


hC|h* 


Vo  + V. 


(2) 
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and 


and 


(3) 


there  the  polarity  Is  defined  by  Harrison 
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&a 


a - 

P 


(Vg  + v„> 


CJJ2 

y 


Here  Vg  and 


are  the  energy  parameters  defined  as: 


< hC(H  | h*  > 

” ! Tc  .a  . 2 

1— |<  h jh  >| 


00 


(5) 


for  the  covalent  energy  where  < hc|  }{|h*  > corresponds  to  hopping  between 
an  anion  hybrid  and  a cation  hybrid  along  the  same  bond. 


[<hC|M|hc>-  <h*|,M|ha>1 
2(1  - j<  h°|  ^*>1^ )' 


for  the  polar  energy  where  < hc| # | hc  > and  < h*|j*(  |h8  > denote  the 

coulomb  integrals  at  the  cation  site  end  the  anion  site  respectively.  In 

8 (c ) 

addition,  there  is  one  more  energy  parameter  , which  is  defined  as 

a matrix  element  of  hamiltonian  between  adjacent  anion  (cation)  hybrids.  This 
la  related  physically  to  the  atomic  s-p  splitting  and  does  not  come  into  play 
In  most  optical  properties. 

The  dielectric  constant,  which  is  the  parameter  of  interest  in  optical 
properties,  is  calculated  by  considering  an  isolated  bond  dipole  with  pn 
applied  electric  field. 


The  field  perturbation  is  treated  by  adding  to  the  polar  energy  the 
potential  ener^  of  the  bond  dipole  in  the  electric  field.  Then  in  the 
ground  state  wave  function  approximation,  calculation  of  the  induced 
polarization  by  a change  of  polar  ity  u ’ through  a change  of  the  polar 

energy  gives  the  dielectric  constant 


K«72(edf  V2 
C°  1 3(V?  + V?)3^ 


(7) 


vbere  5 ■ density  of  electrons,  d • bond  length  and  7 is  a parameter 
defined  in  the  BCM  formulation  relating  microscopic  dipcle  moments  to 
ttacrcscopic  linear  susceptibilities  and  is  taken  to  bt  constant  along  • row 
of  the  periodic  table  except  for  correction  for  carbon  row  atomj. 

20 

To  determine  the  energy  parameters  Harrison  and  Ciraci  appeal  to 
experimental  parameters.  The  principal  optical  absorption  peak  E-.  is 

fc  m\ 

Identified  with  '.he  bonding  and  anti-bonding  gap 


This  identification  has  been  shown  to  be  consistent  with  experiments. 

First,  the  gap  of  the  group  IV  elemental  compounds  gives 

directly  for  each  row  of  the  periodic  table,  since  * 0 in  elemental 
compounds.  The  geometric  mean  is  taken  for  compounds  involving  two  rows. 
Then  using  Eq.  (7)  the  dielectric  constant  gives  7 for  each  row.  A value 
for  7 of  1.08  is  found  for  diamond  and  1.2  for  silicon.  When  row  three 


- m - 


la  reached,  there  are  contributions  from  the  polarizable  d core  electrons 

17 

aa  noted  by  Phillips.  This  is  taken  into  account  by  introducing  a correction 
factor  0 in  addition  to  7 . 7 la  taken  et  the  silicon  value  and  a value  of 

0 ■ l.lfl  is  obtained  for  germanium  and  0 = l.Ul  for  the  tin  row.  With  these 
rallies  of  Vg  and  7 , for  any  binary  compound  is  derived  from  its 

dielectric  constant  using  Eq.  (7)-  This  coopietes  the  summary  of  the  Bond 
Orbital  Model  as  originally  stated  by  Harrison.^ 


m.  CALCULATION  OF  x(2) 

The  second  order  susceptibility  is  defined  by  expanding  the  macroscopic 
polarization  in  a power  series  in  the  field 


_P(E) 


X»>  E 


• : E E + ... 


(8) 


vhere  and  x^2^  *re  the  field  independent  susceptibilities.  In  the 

independent  bond  formulation  these  macroscopic  susceptibilities  are  related 

to  their  corresponding  microscopic  bond  polarizabilities  through  geometrical 

. . 21,22 

factors  as  ' 


V «,t,J  II 


*11 


(9) 


.(2).  i 2^  ,(»)  ,(*>  ,(*>. 

CUK  V s,i,J  aIi  8Ji  aKk  Bijk 


(10) 


vhere  a^  are  the  direction  cosines  between  fhe  macroscopic  crystal  axis  I 
and  the  bond  axis  frame  i . The  index  s identifies  the  different  bonds  in 
the  periodic  lattice  structural  unit  whose  volume  is  V . Tc  apply  the  BCM, 
confine  ourselves  to  sp^-bonded  tetrphedral  compounds  with  the  zinc  blende 
and  wurtzite  structures.  In  the  former  case  with  crystal  cIpss  point  group  “3  m, 
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Eqa.  (9)  and  (10)  reduce  to 


xu)-  ^ (a|l  H 


»A> 


where 


. a 


and  O.  I and  a.  are  the  parallel  and  perpendicular 


* ■ 

coaponents  of  the  bond  polarizability. 

Here  n » 3 JJ/lS  d3  is  the  number  of  tetrahedral  unit  cells  per 
unit  volune.  In  this  esse  x(1)  Is  . sc.lsr,  snd  *£>  U th«  only  non_  ' 

vanishing  tensor  component  in  Voigts  notation.  is  along  the  bond 

axis  (111)  and  with  axially  symmetric  C-y  bonds,  the  independent  non- 
vanishing components  are  o( , - 0^  , ^ - cr^  - 0yy  , and  3 - 3„,, 

" 6xzx  “ syzy  * For  thc  casc  of  wurtzite  structure  the  two  inequivalent 
tetrahedra  can  be  obtained  from  standard  zinc  blende  orientation  by  simple 
rotational  matrices,  the  result  is2^ 


- 


(2)  - - y(2)  . v(2>  2 (2) 

31  y15  *33  * *7=  'it 


It  can  be  shewn  from  an  abinitio  calculation  that  3^  < < 9 ,.2l>  ^ a bond 

chsrec  approach,  3j_  . to  1st  order.  25  ,ceordlRsly>  ue  „egl.ct 
to  in  our  formulation.  ■** 

To  calculate  9 j we  need  to  determine  the  variation  of  the  linear  polar- 
izability a with  an  electric  field.  Physically,  we  expect  a redistribution  of 
charge  density  in  the  presence  of  the  field.  This  manifests  itself  as 
an  induced  dipole  ns  a result  of  -n  additional  charge  transfer.  In  the 

context  of  the  BOM  there  is  a chpnge  in  the  bond  polarity,  a through 

P 
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• change  in  the  bond- polar  energy,  . Harrison  defines  an  effective 
charge  on  an  anion  (in  units  of  -e)  as 

z*  - Ua  - tz  (13) 

p 

vhere  LZ  is  half  of  the  nominal  valence  difference.  Then  the  additional 
charge  incurred  by  the  applied  field  on  a per-bond  basis  is 


6a, 


(14) 


vhere  the  nain  contribution  is  from  the  change  in  polar  energy,  • The 
covalent  energy,  Vg  , is  constant  since  the  bond  length  remains  fairly 
undistorted.  This  is  valid  in  the  optical  transparency  region  above  the  phonon 
lattice  resonance  frequencies,  where  the  lattice  and  the  bonds  stay  relatively 
stationary  compared  to  the  motion  of  the  electronic  charge  cloud.  Here  we 
have  to  stay  below  the  interband  electronic  absorption,  as  the  assumption  of 
a bonding  state  of  the  hybrid  orbitals  implies  that  we  are  dealing  with  valence 
electrons  only.  In  the  transparency  region  the  optical  response  is  fairly 
Insensitive  to  frequency  end  hence  our  dispersionless  approach  should  be  adequate. 

In  the  formulation  of  the  ECM,  when  the  macioscopic  dielectric  constant 
Is  used  to  calculate  the  microscopic  bend  parameters,  the  effects  of  charge 
distribution  scaling  and  the  bond  anisotropy  are  automatically  included  in 
the  parameter  7 . If  we  let  7*  be  the  charge  distribution  scaling  and 
If  ■ a., /a  be  the  anisotropy,  then  as  in  the  derivation  In  BCM,  the  inclusion 

I I • 

of  the  dipole  potential  energy  in  en  electric  field  as  a change  in  the  polar 
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energy  gives  in  terms  of  7* 


*0 


(ed)2 


II  2(V^  + V^)3/2 


(15) 


Bitting  this  into  Eq.  (11)  and  comparing  the  result  with  Eq.  (7) 
we  find  that  72  » 7*2  (1  + 2#f)where  € - 1 + 4 x has  been 

employed.  There  is  no  easy  analytical  calculation  for  the  bond  anisotropy 
By  simple  extrapolation  on  values  obteined  from  Kerr  constants  cn  diatomic 

1 , 26 

molecules,  ve  estimate  K ~ 0.5.  Essentially  the  same  estimate  was 
employed  previously.  With  this  value  for  the  anistropy  factor  and  the 
BOM  value  of  1.2  for  7 , y*'  is  found  to  be  0‘85. 

To  introduce  the  explicit  field  dependence  we  note  that,  the  induced  bond 
dipole  is 


dp  « 7*ed  AZ*  « Oj  | E|  i 


(16) 


As  explained  earlier  an  electric  field  induces  an  additional  charge 
transfer  resulting  in  a change  in  the  polar  energy  . The  second  order 
polarizability  can  be  related  to  its  linear  counterpart  by  QCj  j CE^  ^ ) * 

' ”°f|j+^0f||  + ’**  Now  Aot|  j "P||Ejj  can  be  calculated 
from  Eq.  (15)  through  a change  in  iolar  energy  AV_  . Elimination  of  Av 

3 3 

through  Eq.  (14)  end  Eq.  (16)  in  terms  of  E.i  gives 


'|i  tBiiEn  + 


3 <* 


P. 


ed  (1  - a Y 
P- 


(17) 


To  convert  to  experimental  macroscopi.  quantities,  Eq.  (11)  and 


- 118  - 


Eq.  (12)  give 


y(a).  ki  °p  ,_U).a 

14  «y  (a«  ♦ i)2  (i  - o£) 


eliere  the  tetrahedral  cell  density  (n  « d^)  has  been  used.  Since 

°p  c,n  be  emulated  from  the  linear  dielectric  constant  in  BCM  theory, 

(18)  shows  that  can  be  estimated  essentially  from  Just  the  knowledge 

ot  the  linear  susceptibility  alone.  This  is  an  important  point.  The  independent 
bond  formulation  is  only  approximate  in  the  sense  that  correlation  effects 
between  neighboring  bonds  do  give  rise  to  local  field  corrections  which  are 
not  subject  to  easy  deductive  analysis.  By  appealing  to  the  macroscopic 
linear  susceptibility  as  an  input  parameter,  most  of  these  local  field  effects 
are  accounted  for  in  the  BOM  model. 

To  demonstrate  the  dependence  of  y'  ' ja  polarity  and  to  display 
•one  underlying  physics,  we  transform  Eq.  ( 18)  by  casting  in  terns  of 

polarity  with  the  help  of  Eq.  (7).  Using  the  relation"20  Vg  a l/d2  we  obtain 


y(2)  ~ 

*14  * 


Constant 


a (1  - a2)* 
P P 


where  the  constant  contains  7 , and  is  constant  for  one  row  in  the  periodic 
table. 

• to) 

The  dependence  of  y'  on  polarity  is  shown  in  Fig.  1,' where 
Vd  plotted  against  0^  . The  upper  boundary  corresponds  to  the  fifth 
row,  and  the  lower  one,  the  third  row.  Because  of  differences  in  geometric 


factors,  we  have  only  included  the  case  of  zinc  blende.  Note  the  interesting 


maximum  at  ar^  0.44,  corresponding  to  an  ionicity  ^ « 0.28  in  Phillip's 
model  [f^  = l~(l  * a|,3/2)  . Tang^  and  Levine-**  have  predicted  an  optimum 
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v*lue  of  f^  •»  0.17.  However,  available  experimental  data  seems  to  show  a 

peak  at  around  ■ 0.U5  (f^  ~ 0.3)  in  better  agreement  with  our  results. 

The  physics  of  this  maximum  warrants  scms  comment.  In  the  complete  polarity 

limit  of  0^  - 1,  i.e.  full  ic  -city,  all  excess  valence  electrons  from  the 

cation  are  transferred  to  the  anion.  The  anion  in  turn  binds  these  electrons 

fairly  tightly  in  its  locality,  resulting  in  a big  forbidden  energy  gap 

characteristic  of  ell  ionic  compounds.  Actually,  it  can  be  shown  that  the 

energy  gap  of  tetrahedral  binary  compounds  scales  linearly  with 

This  there  is  no  more  charge  transfer  (dZ  « ) t0  reaP°nd  to  ®ny  applied 

(o\ 

field  and  hence  the  vanishing  of  k . In  the  covalent  limit  of  Of  ■ 0 

P 

ve  approach  the  transition  of  zinc  blende  into  the  elemental  diamond  structure 
and  the  center  of  inversion  symmetry  completely  cancels  all  even  order  of  non- 
linearities.  In  a sense,  the  system  overcompensates  itself  because  the  gain 
in  the  covalent  'bond  charge'  is  completely  cancelled  by  the  perfectly  symmetric 
lattice  environment.  A minor  word  of  caution  stems  to  be  appropriate  here  since 
we  must  remember  that  in  actual  practice,  we  ha\j  to  stay  within  f < 0.78  for 
tetrahedral  coordination.  Any  higher  ionicity  results  in  rocksalt  structure, 
a hexa-coordination  typical  of  most  ionic  compounds. 

To  further  test  the  accuracy  of  our  analysis  we  cnoose  materials  of 
roughly  equal  polarity,  i.e.  dZ  « constant,  and  plot  log^Q  E*gXpt  10'] 
against  log^Q  [d]  in  Fig*  2.  Lines  of  slope  4 are  drawn  and  agreement  is 
good  with  Sq.  (18). 

Thus  for  large  second  order  nonlinearities,  it  is  desirable  to  have  a 
covalency  high  enough  (df^  around  0.44)  to  give  considerable  polarizability 
from  charge  transfer  but  low  enough  such  that  inversion  symmetry  is  not  a 
detriment.  We  also  want  long  bond  lengths,  which,  from  empirical  evidence,  means 
a high  atomic  number.  As  high  covalency  means  proximity  to  Croup  IV,  the  above 
two  criteria  could  be  summarized  In  the  context  of  the  periodic  table  by  saying 
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(9) 

th*t  for  high  y^  we  want  materials  with  smhll  tZ  and  high  atonic  number. 

These  criteria  are  well  borne  out  by  experimental  evidence.  The  III-V 
compounds  are  in  general  acre  nonlinear  than  the  H-V  and  among  the  III-V 
themselves,  the  more  metallic  elements,  InSb  and  InAs,  which  are  from  the 
fifth  rows  and  fifth  and  fourth  rows  respectively,  are  more  nonlinear  than 
OmAs  and  Gap  being  ccnpounds  from  fourth  rows  and  fourth  and  third  rows.  In 
Fig.  3 we  plotted  x^  against  the  average  quantum  number,  n , introduced 

27 

First  by  Mooser  end  Pearson  , verifying  the  above  arguments. 

However , for  any  device  applications  employing  the  second  order  non- 
l^®e*rity,  the  parameter  of  interest  is  the  parametric  gain  coefficient^ 

defined  by 

_*  9.106  I, 

r « r 

8 eJ  n^n^ 

*here  1^  is  the  pump  intensity  and  n^  are  the  refractive  indices  at  oo^  . 

From  the  device  standpoint,  besides  the  magnitude  of  y^  , there  are 
elso  additional  trrnsparency  factors  in  the  frequence  dependence  of  the  gain 
coefficient.  As  one  goes  towards  higher  atomic  numbers  one  loses  ground  in 
the  trensparancy  range  as  determined  by  the  bend  gap  cut-off  frequency. 
Experiment  ally , a material  is  chosen  with  the  smallest  band  gap  consistent 

with  transparency  in  the  wavelength  region  of  interest,  as  coint-d  out  by 

lb 

Icvine.  If  half  the  bandgap  frequency  is  taken  for  both  and  co 
•void  two  photon  absorption,  it  turns  out  that  in  seme  cases  of  high  atomic 
numbers,  e.g.  InAb  and  InAs,  the  increase  in  x^  is  largely  offset  by 
frequency  factors  in  the  gain  expression  due  to  the  narrow  band  gap. 


IP! 


Another  parameter  which  has  proven  to  be  very  useful  as  a scale  for 
nonlinearity  is  the  Miller's  delta  » which  is  defined  by 


2a 

Aijk  *ii 


In  the  context  of  our  model  Miller's  delta  can  be  obtained  by  transforming 
Bq.  (IB)  through  use  of  expressed  in  polarity  a and  bond  length  d . 

We  find  that  A a a d/(l  - a )5^2  . Thus  we  see  that  for  materials  of 

1 jX  p p 

the  same  average  quantum  number  n , where  the  bond  length  is  essentially 
constant  and  for  materials  of  the  same  vertical  family  i.e.  of  constant  polarity 
the  Miller's  delta  should  be  constant.  Essentially  a similar  conclusion  is 
reached  by  Tang  and  Flytzanis,1^  though  they  employed  slightly  different  para- 
meter*, namely  the  normalized  valence  difference  and  the  atomic  radii  difference 
In  Miller's  formulation,  which  is  essentially  an  enharmonic  oscillator  model. 
Miller’s  delta  is  found  to  be  a function  of  the  asymmetric  potential.  Because 
polarity  (charge  transfer)  is  a direct  result  of  the  esymmetry  between  the 
cation  and  anicn,  and  bond  length  also  affects  the  potential,  the  dependence 
of  Miller's  delta  on  the  above  BCM  parameters  is  not  surprising. 

The  main  results  and  parameters  ’.sed  are  summarized  in  Tables  la  and  lb. 

We  have  chosen  the  most  recent  data  or  that  wt 'ch  we  believe  to  be  the  most 
reliable.  Correlation  with  experimental  values  is  in  general  very  good.  It 
•hould  be  observed  that  there  ere  two  general  classes  of  compounds  'here  the 
relative  discrepancy  warrants  some  discussion.  They  8re  the  narrow  band  gap 
materials  for  example  InSb,  InAs  and  GaSb  and  the  copper  compounds. 


1 P2  - 


20 

At  pointed  out  in  the  original  formulation  of  BCM,  the  use  of  the 

peak  for  the  bonding  and  anti  bonding  gap  works  especially  well  for 
the  large  band  gap  materials.  For  narrow  band  gap  materials,  the  minimum  optical 
g*p  could  be  almost  an  order  of  magnitude  smaller  than  the  gap.  For 

example  in  InAs  EQ  » 0. 33  and  ■ 3*9  eV  so  one  might  expect  the  r 

point  in  the  Brillouin  2one  to  modify  the  BCM  energy  parameters.  As  a result, 
the  application  of  the  BGM  to  these  compounds  is  not  as  successful. 

For  the  copper  compounds,  it  is  known  that  the  d-bands  art*  strongly 
hybridized  with  the  valence  s and  p bends.  These  delocalized  d electrons 
at  these  noble-metal  cations  have  been  found  experimentally  to  give  an 
anomalous  contribution  to  the  second  order  nonlinearity. effect 

^Ji  a / 

la  discussed  by  Levine  in  his  bond  charge  model  and  by  Tang  in  terms  of 
an  charge  redistribution  and  a reversal  of  the  direction  of  the  bond  dipole. 
8ince  the  BCM  deals  with  the  sp^  electrons,  the  d-electrons  are  accounted 
for  phenomenologically  and  the  formulation  overestimates  the  nonlinearity  of 
these  compounds.  If  we  stay  away  from  these  noble  metal  compounds  and  consider 
compounds  where  there  is  no  appreciable  hybridization  of  the  d-bands,  then 
the  BOM  predictions  agree  well  with  experiment. 

IV.  SUMMARY  AND  COMMENTS 

Using  a more  predictive  theory,  the  BCM,  we  have  calculated  the  second 
order  susceptibility  for  tetrahedral  compounds.  Agreement  with  experiment 
Is  good,  especially  in  light  of  the  simplicity  of  the  model. 
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The  approach  is  self-consistent  as  the  linear  susceptibilit/  is  used  to 
calculate  the  second  order  susceptibility.  Compared  to  other  phenomenological 
models,  our  approach  is  more  self-contained,  as  no  adjustable  parameter  is  used 
and  all  parameters  are  defined  in  the  BOM. 

The  physical  origin  of  the  nonlinearity  is  seen  to  arise  from  the 
charge  transfer  as  a result  of  the  asymmetric  pciar  energy  between  the 
anion  and  the  caticn.  Chemical  trend  considerations  reveals  an  optimum 
value  of  polarity  for  a maximum  second  order  nonlinearity  in  zinc  blende  and 
wurtzite  structures.  This,  together  with  the  bond  length  dependence  gives 
a rule  of  thumb  estimate  on  the  choice  of  elements  in  the  periodic  table 
which  will  result  in  tetrahedral  compounds  of  high  second  order  nonlinearity. 

It  is  found  that  for  these  tetrahedral  structures,  it  is  desirable  to  have 
high  covalency  end  a high  metallicity  consistent  with  transparency  region  of 
interest.  This  means  that  one  should  stay  close  to  Group  IV  elements  and 
high  atonic  numbers.  We  include  the  parametric  gain  paremeter  in  our  discussion 
showing  that  the  band  gap  is  again  important  through  the  short  wavelength 
frequency  cut-off  in  the  transparency  range.  The  bond  orbital  model  is  also 
quantitatively  and  qualitatively  related  to  the  Miller's  rule  through  the 
concept  of  polarity. 
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TABLE  I* 


COMPARISON  BETWEEN  THEORETICAL  Y^  (CALC.)  AND  EXEERIMEZITAL  (EXPi.  ) 

TOGETHER  WITH  THE  RELEVANT  ECM  PARAMETFBS  "Jl  ZINC  BLENDE  STRUCTURES 


TABLE  lb 


Compound  Bond  length 
d [ A ] 


X^(c*!)  *^(exPfc) 


1.88 

0.39 

0.45 

0.93 

1.98 

O.70 

0.24 

0.88 

2.34 

0.73 

0.33 

2.64 

2.53 

0.77 

0.33 

3.7k 

2.63 

0.77 

0.38 

5*35 

0.78  h 


0.33  i 


1.8  .f 


3.8  J 
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FIGURE  CAPTIONS 


(2)  U 

1*  The  scaling  of  /d  with  °P  t polarity.  The  upper  boundary 

corresponds  to  aaterials  from  the  fifth  row  and  the  lower  one  from  the 
third  row. 
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The  dependence  of 
constant  polarity. 


on  bond  length  d 


for  materials  of  roughly 
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A plot  of  XjY  versus  n , the  average  quantum  number,  showing 
tbe  dependence  within  a vertical  family  in  the  periodic  table. 
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